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Abstract. 

The measurements of reactions n~p^ n~n~^n at 17.2 GeV/c and n~^n^ 
7T^'K~p at 5.98 and 11.85 GeV/c made at CERN with polarized targets provide 
a model-independent and solution-independent evidence for a narrow scalar state 
cr(750). The original minimization method and the recent Monte Carlo method 
for amplitude analysis of data at 17.2 GeV/c are in excellent agreement. Both 
methods find that the mass distribution of the measured amplitude |<S'pS with re- 
coil trans versity "up" resonates near 750 MeV while the amplitude \S\'^T, with recoil 
transversity "down" is large and nonresonating. The amplitude \S\'^T, contributes 
as a strong background to S-wave intensity Is = {\S\'^ + \S\^)Ti and distorts the 
determinations of a resonance parameters from Is- To avoid this problem we per- 
form a series of Breit-Wigner fits directly to the measured distribution |S'pE. The 
inclusion of various backgrounds causes the width of (t(750) to become very nar- 
row. Our best fit with t-averaged coherent background yields = 753 ± 19 MeV 
and Tcr = 108 ± 53 MeV. These values are in excellent agreement with EUis-Lanik 
theorem for the width of scalar gluonium. The gluonium interpretation of a"(750) 
is also supported by the absence of a (750) in reactions 77 tttt. We also show 
how data on polarized target invalidate essential assumptions of past determina- 
tions of TTTT phase shifts which explains the absence of cr( 750) in the conventional 
phase shift Sq. We examine the interference of (t(750) with /o(980) and find it has 
only a very small effect on the determination of a (750) mass and width. The data 
on amplitude I^SpE in the mass range of 1120-1520 MeV show existence of a scalar 
resonance /o(1300) with a mass of 1280 ± 12 MeV and a width of 192 ± 26 MeV. 
We point out that the study of production processes on the level of spin amplitudes 
measured in experiments with polarized targets may reveal new hadron structures 
and open new physics beyond the standard QCD quark model. 
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I. Introduction 

Science is an ongoing interplay between our ideas and our experiences in the 
real Universe. We cannot expect to make a progress in our understanding unless 
we can show that what we thought we knew is in some sense incomplete or even 
wrong. ^ Amplitude analyses of hadronic reactions provide a clear confirmation of 
this logic of scientific discovery. 

In 1972, L. van Rossum and his group at Saclay measured recoil nucleon po- 
larization in TT^p TT^p reactions^ at CERN-PS. Their results closed the set of 
complete measurements of ttN — > ttN reactions at 6 GeV/c and allowed them to 
make the first model independent determination of hadronic amplitudes directly 
from scattering data.^ Failures of Regge models to correctly predict polarization 
and other spin observables were traced to the wrong structure of their amplitudes. 
Some Regge models were revised and constrained to reproduce the experimental 
ttN amplitudes at 6 GeV/c. 

In 1978 the Saclay group reported^ the first amplitude analysis of KN and 
KN charge exchange reactions, also at 6 GeV/c. The structure of ^2-exchange 
amplitudes was found different from all revised Regge models. A more difficult 
amplitude analysis of pp — > pp at 6 GeV/c was made possible by polarized proton 
beam at Argonne ZGS.^ The results reported^ in 1985 by A. Yokosawa and his 
group at Argonne confirmed our lack of understanding of hadronic reactions at the 
level of amplitudes. So far, no new revisions of Regge models were even attempted. 
When the first measurement of a complete set of observables in elastic scattering of 
protons on ^^C carbon nucleus was reported^ in 1981, it invalidated the standard 
nonrelativistic analysis in favour of relativistic approach to nuclear physics.^ 

The work by Saclay, Argonne and other groups firmly established the need 
for experimental knowledge of hadronic amplitudes. Efficient acquisition of this 
knowledge in two-body reactions has been hampered by the difficulty of measuring 
the recoil nucleon polarization. Dispersion relations were used in an effort to obtain 
hadronic amplitudes from incomplete data.^ However, the situation is different in 
pion production reactions ttN — > ttttN and KN — > KttN. In 1978, Lutz and Rybicki 
showed^° that measurements of pion production in meson-nucleon scattering on 
transversely polarized target yield in a single experiment enough observables that 
almost complete and model independent amplitude analysis can be performed. 

Amplitude analyses of pion production reaction such as nN —>■ ttttN or 
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NN — > ttNN are important for two special reasons. First, these reactions pro- 
vide information about unnatural exchange amplitudes which are not accessible in 
two-body reactions. Second, such amplitude analyses enable us to study resonance 
production on the level of spin-dependent amplitudes rather than spin-averaged 
cross-section d? a/dmdt. 

The high statistics measurement of 7r~p n~'K~^n at 17.2 GeV/c at CERN- 
PS on unpolarized target was later repeated with a transversely polarized pro- 
ton target at the same cncrgy^^ Model independent amplitude analyses were 
performed for various intervals^^ of dimeson mass of small momentum trans- 
fers —t = 0.005-0.2 (GeV/c)^ and over a large interval^^ of momentum transfer 
-t = 0.2- 1.0 (GeV/c).2 

Additional information was provided by the first measurement of Tr'^n — > 
TT+TT^p and K'^n — > K'^7r~p reactions"*^^'-^^ on polarized deuteron target at 5.98 and 
11.85 GeV/c, also done at CERN-PS. The data allowed to study the t-evolution of 
mass dependence of moduli of amplitudes. -"^^ Detailed amplitude analyses^^'^-*^ deter- 
mined the mass dependence of amplitudes at larger momentum transfers — t = 0.2- 
0.4 (GeV/c).2 

Amplitude analyses of 7r~p — > TT~TT'^n and ir'^n — > tt'^ti~p reactions were re- 
cently repeated^^'^^ with special attention paid to error propagation and selection 
of physical solutions. This work was motivated by the emerging evidence from 
previous analyses for a new scalar resonance (t(750). 

All amplitude analyses^^"^^ of pion production on polarized targets found a 
clear evidence for large and nontrivial unnatural Ai-exchange amplitudes in the di- 
pion mass range from 400 to 1800 MeV. This experimental finding is very important 
since previously the Ai exchange amplitudes were assumed absent. In particular, 
all determinations of tttt phase shifts from unpolarized data on tt~p — > Tr^Tr+n 
are based on the assumption of vanishing Ai exchange amplitudes. Without 
this assumption the determination of tttt phase shifts cannot even proceed. The 
existence of large and nontrivial A\ exchange amplitudes in ttN — > tt'^tt~ N reac- 
tions casts a serious doubt about the reliability of the conventional tttt phase shifts. 
The assumption of absence of Ai exchange amplitudes means that pion production 
in ttN — > 'K'^'K~N does not depend on nucleon spin. What the measurements of 
ttN TT'^TT~N on polarized targets found is that pion production depends strongly 
on nucleon spin. The dynamics of pion production is therefore not as simple as was 
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assumed in the past determinations of tttt phase shifts. 

The existence of exchange is also crucial for our understanding of spin 
structure of the nucleon. The measurements of the cross-section asymmetry using 
longitudinally polarized lepton beams on a longitudinally polarized nucleon targets 
determine nucleon spin dependent structure function gi{x, Q'^). The measurements 
of gi on proton and neutron targets at CERN and SLAC has fascinating implications 
for the internal structure of the nucleon. These analyses of nucleon spin structure 
depend on the behaviour of gi for a; ^ which is controlled by Ai exchange (see 
eq. (4.2.23) of Ref . 32). 

Another important finding of measurements of 7r~p Tr^Tr+n and Tr+n 
TT~^7T~p reactions on polarized targets is the evidence for a narrow scalar state / = 
0++(750). 

The first published evidence for this state in 7r~p — > tt^tt+tt, is presented in 
Fig. 2 of Ref. 14 from 1979. The CERN-Munich data on ;S-wave intensity Is show 
a clear bump around 750 or 800 MeV (depending on a solution). The CERN-Munich 
group refers to this bump as a resonance e(800) but makes no Breit-Wigner fits. In 
1979, Donohue and Leroyer published^^ an analysis of the CERN-Munich data on 
polarized target and made the first claim that the data show existence of a narrow 
resonance which they called e(750). 

Our Saclay group measured yr+n — > tt+tt"!) on polarized deuteron target at 
5.98 and 11.85 GeV/c at CERN-PS. The measured 5"- wave intensity Is showed a 
narrow resonant structure at 750 MeV at both energies (see Fig. 10 of Ref. 20). To 
verify our results on Is in 7T~^n tt^tt~p reaction, we analysed the CERN-Munich 
data at 17.2 GeV/c using our computer programs. All four solutions for 5'-wave 
intensity Is at all 3 energies showed a narrow resonant structure around 750 MeV. 
We reported these results in 1992 in Ref. 22 and made the claim for existence of a 
narrow scalar state (7(750). However, later a numerical error was found in the final 
step of calculation of S-wave intensity Is at 17.2 GeV/c. When corrected, only 
2 solutions for Is at 17.2 GeV/s show a clear narrow structure around 750 MeV 
while the other two solutions have a broader behaviour. Thus there appeared some 
difference between the results at low momentum transfer data at 17.2 GeV/c and 
higher momentum transfer data at 5.98 and 11.85 GeV/c. 

In Ref. 22 we presented analytical solutions of our amplitude analyses which 
included the unphysical solutions in many (m, t) bins. To deal with the unphysical 
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solutions and to improve our error analysis we used a Monte Carlo method of 
amplitude analysis in Ref. 23. A clear signal for a narrow (t(750) emerged from this 
improved analysis. This evidence for a narrow cr(750) state was recently confirmed 
in a new measurement of 7r~p — > Tr^Tr"*"??- on polarized target at 1.78 GeV/c made 
at ITEP accelerator in Moscow and reported^^ in 1995. 

To describe our new results, compare them to the CERN-Munich analyses, 
and to outline the program of this paper, we first need to introduce our notation 
and discuss the two accepted methods for amplitude analysis of data on polarized 
targets. 

For masses below 1 GcV the dimeson system is produced predominantly in 
spin states J = (5'-wave) and J — 1 (P-wave). The experiments yield 15 spin- 
density-matrix (SDM) elements describing the dimeson angular distribution. These 
observables and the cross-section dPa/dmdt = S can be expressed in terms of two 
5- wave and six P-wave nucleon transversity amplitudes. -"^^'^^ In our notation, the 
two 5'-wave amplitudes are S and S. The six P-wave amplitudes are L, L, t/, U 
and Ar,iV. The amplitudes 'A = S,L,U,N and A = S,L,U,N correspond to 
recoil nucleon transversity "up" and "down" relative to the scattering plane. The 
amplitude analysis works with normalized amplitudes using a normalization 

l-Sp + l^p + |Lp + |Lp + + |t7p + |Ar|2 + |iV|2 = 1 (1.1) 

The unnormalized amplitudes then are \A\'^T, and We also define partial- 

wave intensity 

Ia = {\A\^ + \A\^)^ (1.2) 

where A = S,L, U, N. There are two solutions for the moduli and independent 
two solutions for the moduli \A\'^. Hence, there are 4 independent solutions for the 
partial wave intensities I a which we label lA{i, j = 1,2 with indices i and j 
referring to the two solutions for \A\^ and \A\'^, respectively. 

Amplitude analysis expresses analytically^°'^° the eight normalized moduli and 
the six cosines of relative phases of nucleon transversity amplitudes in terms of mea- 
sured SDM elements. There are two similar solutions in each (m, t) bin. However, 
in many (m, t) bins the solutions are unphysical: either a cosine has magnitude 
larger than one or the two solutions for moduli are complex conjugate with a small 
imaginary part. Unphysical solutions also complicate the error analysis. 
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The occurrence of unphysical solutions is a common difficulty in all amplitude 
analyses. Two methods are used to find physical solutions and determine their 
errors. They are (a) minimization method and (b) Monte Carlo method. 

In the method one minimizes a function 

2 Obsi(meas) - Obsi(calc) 2 o\ 

X =1^[ X, ] (^-^^ 

where Obsj (meas) are the experimentally measured quantities, Aj are their ex- 
perimental errors and Obsj (calc) are corresponding expressions in terms of the 
amplitudes (moduli and cosines of relative phases) . The analytical solutions for the 
moduli and cosines serve as initial values. This x^ method was used in all CERN- 
Munich analyses^^"^^ of 7r~p — > tt^tt+p at 17.2 GeV/c. Since the two analytical 
solutions (initial values) are very close, the x^ method leads to a unique solution in 
many (m, t) bins. A particular exception is the mass range below 900 MeV. More 
recently the x^ method was used in direct reconstruction of amplitudes of pp elastic 
amplitudes from 0.8 to 2.7 GeV using polarized data obtained at SATURN II at 
Saclay.^^ 

The basic idea of Monte Carlo method is to vary randomly the input SDM 
elements within their experimental errors and perform amplitude analysis for each 
new set of the input SDM elements. The resulting moduli and cosines of relative 
angles are retained only when all of them have physical values in both analytical 
solutions. Unphysical solutions are rejected. The distributions of accepted moduli 
and cosines define the range of their physical values and their average value in 
each (m, t) bin. The Monte Carlo amplitude analysis of Ref. 23 is based on 30,000 
random variations of the input SDM elements. The Monte Carlo method was first 
used in 1977 in an amplitude analysis^^ of pp elastic scattering at 6 GeV/c and 
later in an amplitude analysis^^ of reactions 7r~p K'^K~n and 'K~p KgKgn 
at 63 GeV/c. In his review paper,^^ F. James advocates the use of the Monte Carlo 
method as perhaps the only way to calculate the errors in the case of nonlinear 
functions which produce non-Gaussian distributions. The method has the added 
advantage that it can separate the physical and unphysical solutions and that it 
can retain the identity of the two analytical solutions. 

The results for the two solutions for the unnormalized moduli of S-wave am- 
plitudes l/SpS and obtained by Monte Carlo amplitude analysis of CERN- 
Munich data at 17.2 GeV/c are shown in Fig. 1. We find that both solutions for the 



7 



amplitude \S\^T, resonate around 750 MeV while both solutions for the amplitude 
jS'pS show non-resonant behaviour and increase with dipion mass m. 

The results for jS'pS and j/SpS obtained by minimization method using 
the same CERN-Munich data at 17.2 GeV are shown in Fig. 2. We again find that 
both solutions for the amplitude resonate around 750 MeV and that both 

solutions for the amplitude |S'pE show non-resonant behaviour and increase with 
dipion mass m. The comparison of Fig. 1 and Fig. 2 shows that Monte Carlo method 
and minimization method are also in excellent numerical agreement. However, 
the amplitudes obtained by Monte Carlo method show a considerably smoother 
behaviour which, as we shall see later, gives much lower values in Breit-Wigner 
fits. 

The unnormalized moduli \S\'^'Z and |5'pE in Fig. 1 and 2 were calculated using 
S = (fa/dmdt from Fig. 12 of Ref. 11. 

At this point we note that the Fig. 2 is based on Fig. 10 of Ref. 13 and Fig. VI- 
21 of Ref. 12 (to resolve error bars). The authors of these papers present only 
normalized moduli and and consequently did not see the resonant behaviour 
of unnormalized amplitude jS'pS. The resonant behaviour of amplitude j/SpS at 
750 MeV went also unobserved in the subsequent analysis in Ref. 14 which was 
using polarized data in 40 MeV bins in the mass range from 600 to 1800 MeV. It is 
possible to reconstruct the amplitudes iS'pS and \S\'^Tj from the information given 
in Ref. 14. As we shall see in Section V (Fig. 8) both solutions for |5'pE resonate 
below 900 MeV while both solutions for |5'pE are nonresonating, in agreement with 
Figures 1 and 2. It is interesting to note that the evidence for a narrow resonance 
a"(750) was hidden in the very first analyses of CERN-Munich data (Ref. 12, 13 
and 14) and was recognized by the present author some 16 years later in connection 
with the present work. 

The aim of the present work is a more reliable determination of mass and 
width of cr(750) resonance from model independent amplitude analyses of CERN- 
Munich data on 'K~p — > 'K~'K'^n on polarized target at 17.2 GeV/c. There are three 
important issues that we address in the process. 

The first issue in the question which mass distribution should be used for Breit- 
Wigner fits to determine the resonance parameters of a-(750) state. The previous 
CERN-Munich analyses^'^"^^ fitted a Breit-Wigner formula to partial wave intensi- 
ties, and we followed the same procedure in Ref. 23. However, the S'-wave intensity 
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at lower momentum transfers at 17.2 GeV/c shows a clear resonant structure only 
in solutions 15(1, 1) and 15(2, 1) while the solutions Is{^, 2) and 15(2, 2) lack suffi- 
cient decreases of Is above 800 MeV to indicate a narrow resonance. This behaviour 
in Is = (I'S'p + |'S'p)S is caused by the large and nonresonating amplitude |<S'pS. 
The amplitude \S\'^T, thus behaves as a large and nonresonating background to 
the resonating amplitude \S\^T, and this distorts the determination of resonance 
parameters of a (750) from Breit-Wigner fits to Is- To avoid this problem, it is 
necessary to perform Breit-Wigner fits directly to the resonant mass distributions 
|S'pS. Both solutions to |5'pS resonate and the evidence for cr(750) is thus entirely 
solution independent. 

The second issue is which resonance shape formula is to be used in Breit- 
Wigner fits to l^pE. The previous analyses ^^'^^ used the Pisut-Roos shape 
formula which multiplies the standard Breit-Wigner formula (with a phase space) 
by an additional mass dependent factor F = (2 J + l){m/q)'^. In their analysis of 
ttN — > ttttN reaction amplitudes, Pisut and Roos assumed absence of Ai exchange 
amplitudes and assumed that the mass dependence of pion production amplitudes 
is given by tttt scattering amplitudes. The partial wave expansion of tttt amplitudes 
then directly leads to the additional factor F. However, because of the existence 
of large and nontrivial Ai exchange amplitudes in ttN — > tt'^tt~N reactions and 
because there is no proof that the mass dependence of ttN — > titiN production 
amplitudes is really described by the mass dependence of partial waves in tttt scat- 
tering, it is useful to perform Breit-Wigner fits to |5'pE mass distribution using 
the standard'^'^ phenomcnological shape formula in which the Pisiit-Roos factor is 
absent (i.e. -F = 1) to see if there are differences in the determination of resonance 
parameters of (t(750). The comparison of fits using both shape formulas finds only 
small differences. 

The third issue is the question of background in the resonant mass distribution 
Nonresonating background comes e.g. from the isospin 1 = 2 contribution 
to l^pS. While it is difficult to exactly parametrize the unknown background, we 
estimated the background contribution using 3 different approaches. In each case 
we find that inclusion of background leads to a significant reduction of the width 
of (t(750) to somewhere around 100 MeV. Background is obviously important for 
the width determination of a (750) and thus to our understanding of the constituent 
structure of the a (750) state. We also examine the interference of ex (750) with 
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/o(980) and find it has only a small effect on the mass and width of cr(750). 

The paper is organized as follows. In Section II we review the basic formalism. 
In Section III we derive the Pisut-Roos shape formula and describe the phenomeno- 
logical shape formula for Breit-Wigner fits. In Section IV we present our fits to the 
measured resonating amplitude I^SpS and describe our approaches to inclusion of 
coherent background. In Section V we study the interference of cr(750) with /o(750) 
and also perform fits in the broad mass range 600-1520 MeV which show evidence 
for a scalar resonance /o(1300). In Section VI we present our fits to 5'-wave intensity 
Is in n~p 'K~'K^n at 17.2 GeV/c and in tt+ti n~^'n:~p at 5.98 and 11.85 GeV/c. 
In Section VII we review the assumptions of past determinations of nn phase shifts 
and show how they are invalidated by the data on polarized targets. This explains 
the absence of narrow cr(750) state in conventional phase shift Sq. In Section VIII 
we answer critical questions concerning the evidence for a narrow a (750). In Section 
IX we propose to identify the cr(750) state with lowest mass gluonium 0'^'^{gg) and 
discuss theoretical and experimental support for this interpretation of cr(750). In 
Section X we comment on the significance of studying hadron production on the 
level of spin amplituds and suggest that such studies may reveal new physics beyond 
the conventional QCD quark model of hadrons. The paper closes with a summary 
in Section XI. 
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II. Basic Formalism 

A. Phase space and amplitudes . 

Various aspects of phase space, kinematics and amphtudes in pion production 
in ttN ttttN reactions are described in several books. ^-"^"^^ In our discussion we 
will follow often the book by Pilkuhn.^^ 

Consider reaction a+b — > 1 + 2+3 such as 7r~p — > tt^tt+tt- with four-momentum 
conservation 

P =Pa+Pb=Pl+P2+PS=Pd+P3 (2.1) 

where Pd = Pi + P2 is the dipion momentum. The spin averaged cross-section is 
given by 

= Flu!cM(2.„ + l)k+l) ,^ |M..,„.J^dLip.3 (2.2) 

where the flux 



Flux(s) = 4.^{pa-pby-mlml = 4MP^iab (2.3) 

with rua = A* mass of pion and 777,5 = M mass of proton. The Xp and A„ are the 
proton and neutron helicities. The dipion state does not have a definite spin and 
helicity. The Lorentz invariant phase space is defined as 

dLips3 = i2nr6HP -p,-p2- Ps) n ^2nn^E,) ^^"^^ 

We will work with the usual kinematic variables of cm. energy squared s, momen- 
tum transfer t, dipion mass m and angles 9, (f) describing the angular distribution 

of n~ in the n~n~^ rest frame. Hence 

s= {Pa+ Pbf = {Pi +P2+ Psf 

t= {Pa- PdY = {Pa -Pl- P2f (2.5) 

m'^ =pI = {pi+P2f 

Following the procedure described on pp. 18-19 of Ref. 40 and using dm^ — 2mdm, 
we can write 

dLips^{P,pi,P2,P3) = q{m^)G{s)dmdtdil 
where q is the pion momentum in the c.m.s. of the dipion system 



q(rr?) = ^0,25m^ - li'^ = -Jl-A'' (2.6a) 

2 V m 
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and the energy dependent part of phase space 



where 



Hence 



where 



G(s) = — ^ — , ^ (2.6b) 



A(s, M^) = [s-(ii + Mf] [s-ifi- Mf] 

5;^ = ^^EiM^..oA.(M.™.M)i^ (2.7) 



Flux(s) 

To obtain dipion states of definite spin J and hehcity A, we expand 

00 + J 

Mx„,ox, = Y^ /(2J+l)M/,„,OA,(«,^,H^^^oWe^^'^ (2-8) 

J=0 A=-J 

We now integrate \M\^^q\^\'^ over dfl. Using orthogonality relations for the d- 
functions and spherical harmonics, we obtain for the reaction cross-section 

^ = g(m^)i^(.) X: E \ML.,oxAsArn)f (2.9) 

t7=0 A, , Ap 

We now define 

^AA„,OA, = V^y/K(s)M^x^,ox, (2-10) 

We will also consider only the J = (S'-wave) and J = 1 (P-wave) contributions 
for dipion masses m below 1000 MeV. With a notation 

S = d^a/dmdt (2.11) 

we now define normalized helicity amplitudes with definite t-channel naturality 

Hi,^^,+ = ^iNo±Uo)V^ (2.12) 
12 



In (2.12) n= \Xp — A„| = 0, 1 is the nucleon helicity flip. At large s, the unnatural 
helicity nonflip amplitudes Sq, Lq,Uo and the unnatural helicity flip amplitudes 
Si, Li, Ui exchange tt and Ai quantum numbers in the t-channel, respectively. Both 
natural exchange amplitudes A^o and Ni exchange A2 at large s. 

The amplitude analysis of data on polarized targets is performed using nor- 
malized recoil nucleon transversity amplitudes defined as 

S=^iSo + iSi) s=^{So- iSi) 

1 -1 
L = -j={Lo + iLi) L^—{Lo-iLi) 

U=^{Uo + tUi) Tl=^{Uo-iUi) (2.13) 

N = ^{No- iNi) N=l=(No + iNi) 

The amplitudes S, L, U, N and S, L, U, N correspond to recoil nucleon transver- 
sity "down" and "up", respectively.^®'^° The "up" direction is the direction of nor- 
mal to the scattering plane defined according to Basel convention by x^tttt where 
Ptt and Ptttt are the incident pion and dimeson momenta in the target nucleon rest 
frame. 

The normalized amplitudes satisfy conditions 



J2\Sn\^ + \Ln\''+\Un\'' + \Nn\'' = l 



n=0 

l-Sp + \S\'^ + |Lp + |L|2 + \Uf + \Uf + \Nf + \Nf = 1 (2.14) 
We now define spin-averaged partial wave intensities for amplitudes A = S, L,U, N 

Ij, = (|Ao|' + |Ai|2)S = (|A|2 + |Z|2)S (2.15) 

Obviously 

^ = -^ = Is + Il + Iu + In (2.16) 
dmdt 

^Prom the point of view of Breit-Wigner fits to various mass distributions in 7r~p 
'K~'K'^n (S,/yi, etc.) the important point is that the part of the phase space 

which depends on the dipion mass m is simply the cm. pion momentum q given 
by (2.6a). 

13 



B. Spin observables and amplitude analysis. 

For invariant masses below 1000 MeV, the dipion system in reactions tiN 
Ti^T{~N is produced predominantly in spin states J = (5'-wave) and J = 1 
(P-wave). The experiments on transversely polarized targets then yield 15 spin- 
density-matrix (SDM) elements describing the dipion angular distribution. The 
measured SDM elements are^^'^^ 

Pss + Poo + 2pii, Poo - Pii, Pi-i, (2.17a) 

Repio, Repis, Repos, 

pL + P^o + 2p!i,P^o-P^i,pLi, (2.176) 

Imp^ _ 1 , Imp^o , Ixnpl^ . {2.11 c) 

The SDM elements (2.17a) are also measured in experiments on unpolarized targets. 
The observables (2.17b) and (2.17c) are determined by the transverse component 
of target polarization perpendicular and parallel to the scattering plane ttN — > 
{7r'^7r~)N, respectively. The SDM elements (2.17) depend on s,t, and m. There 
are two linear relations among the matrix elements in (2.17): 

Pss + Poo + 2pii = 1, 

pL + P^o + 2p5'i = A (2.18) 

where A is the polarized target asymmetry. 

The data analysis is carried out in the t-channel helicity frame for the tt+tt" 
dimeson system. The helicities of the initial and final nucleons are always defined 
in the s-channel helicity frame. Using the notation of (2.17) and (2.18), the first 
group of equations represents the sum of SDM elements (2.17a) and (2.17b): 

ai = l[l + A] = |5p + |L|2 + + |]v|2, 

ai = [(poo - Pii) + (P^o - Pn)] = 2|L|' - \U\^ - 

as = [Pi-i + Pi-i] = \W- \U\^ (2.19a) 

^[Repio + Rep 
v2 



a4 



05 = 





\L\ 


COs(7l;7), 


= \u\ 


\S\ 


cos(7s[/), 



ae = ^[Repos + Rep^.] = \L\\S\ cos{^sl). (2.196) 
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Similar equations relate the difference of SDM elements to amplitudes of opposite 
transversity. The second group of observables is defined as 

a^ = \[l-A] = \S\'' + |L|2 + \U\^ + \N\\ 

a2 = [(poo - Pii) - {plo - Pn)] ^m^-\W- 

03 = [pi-i - p\-i] = |iV|' - \U\\ (2.20a) 

^4 = -^[Repio - Repfo] = \U\\L\ cos{^lu)^ 

_ 1 

05 = -^[Repis - RepfJ = \U\\S\ cosijsu)^ 

ae = ^[Repos - Rep[j,] = |Xp| 003(75^), (2.206) 

In Eqs. (2.19b) and (2.20b) we have introduced explicitly the cosines of relative 
phases between the nucleon transversity amplitudes. 

The SDM elements (2.17c) form the third group of observables/°'^° which is 
not used in the present amplitude analysis. 

The first group (2.19) involves four moduli jS'p, |Lp, and |A^p and three 
cosines of relative phases cos(75i), cos(75;7), and cos(7^[/). The second group (2.20) 
involves the same amplitudes, but with opposite nucleon transversity. Analytical 
solution for these amplitudes in terms of observables was derived in Ref. 10 and 20. 
For the first group one obtains a cubic equation for |Lp = 

ax^ + bx^ + cx + d=0, (2.21) 

with coefficients a, 6, c, d expressed in terms of observables ai,i = 1, 2, . . . , 6. The 
remaining moduli and the cosines are given by the expressions 



2 

1 



l^"!^ = (ai + 02) — 3|L 
|C/|2 = |L|2-l(a2+a3), 

\W = \Lf-l{a2-as), 

, , a4 

COS(7Li7) 



L\\U\ 



cosi'y su) = 7^r|n ' cos (751.) = T^TT^- (2-22) 



The solution for the second group (2.20) is similar. 
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The analytical solutions of the cubic equation (2.21) are given in Table I of 
Ref. 20. One solution of (2.21) is always negative and it is rejected. The other 
two solutions are generally positive and close. However, in a number of (m, t) bins 
we get unphysical values for some cosines and in some cases also negative moduli 
of amplitudes. In some (m, t) bins the mean values of input SDM elements yield 
complex solutions for |L|^ or |Lp or both (with positive real parts). To filter out 
the unwanted unphysical solutions and to determine the errors on the amplitudes 
and their average values, one can use either the minimization method or the 
Monte Carlo method. The results of amplitude analyses of 7r~p n~'K~^n at 17.2 
GeV/c for dipion masses in the range 600-900 MeV are given in Ref. 13 for the 
minimization method and in Ref. 23 for the Monte Carlo method. The two methods 
are in excellent agreement. In Ref. 23 we also present the Monte Carlo amplitude 
analysis of the reaction Tr+n tt^'k~p at 5.98 and 11.85 GeV/c using the Saclay 
data^^ at larger momentum transfers —t = 0.2 — 0.4 (GeV/c)^ and dipion masses 
in the range 360-1040 MeV. 

III. Resonance shape formulas. 

A. Pisiit-Roos shape formula. 

Before we review the Pisut-Roos derivation of their resonance shape formula, 
we first recall some properties of partial waves in elastic scattering of scalar particles. 
The T-matrix amplitude of isospin I has partial wave expansion 

oo 

T^(s,cos6') = 87r^(2L + l)T/(s)PL(cos6') (3.1) 

The unitarity in elastic scattering then requires (see Ref. 40, pp. 38-40) that T[ 
has a form 

T[{s) = ^ su.Sle< = ^ J, . (3.2) 

q q cotg()[ - I 

where q is the pion c.m.s. momentum and 5£ is the corresponding phase shift. 
Notice that the factor \^/q is induced by the unitarity alone. At a resonance rriR, 
the relativistic Breit-Wigner formula for then reads 

T[ = —-, . , 3.3 

^ q {s-ml)+imRr{s) ^ ' 

where r(s) is an energy dependent width. 
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Let us now return to pion production process 7r~p — > 7r~7r'''n and to amplitudes 
^xXn oXpi^i ^)"^) defined in (2.8). In their analysis,^^ Pisut and Roos assumed that 
the following amplitudes vanish for all J: 

Mo^+,o+ = (3.4a) 

M^i+,o+ = = (3.46) 

The conditions (3.4) mean that all Ai-exchange amplitudes vanish and that the 
natural yl2-exchange amplitudes also vanish. Only pion exchange amplitude Mq_^ q_ 
contribute and they have a general form 



M^+^o-{s,t,m) = Q{s,t)y/2J +lT\m)^Jf{m)+Mi{s,t,m) (3.5) 
where T"^ {m) are the tttt tttt partial wave amplitudes with isospin decomposition 

T"^ = T/^i for J odd (3.6) 

2 1 

J.J ^ _yj ^ -TLo for J even 

3 " 3 

in reaction 7r"'"7r~ — > 7r"''7r~. In (3.5) the function /(m) is a phenomenological 
function that is supposed to account for absorption, and initial-state and final-state 
interactions. In practice one puts /(m) = 1. The function Q{s,t) factorizes the s- 
and t-dependence. The term M^{s,t,m) is a background. 

Taking into account the factor qK{s) in (2.9) and the equation (3.3), Pisut and 
Roos arrive at a resonant parametrization of reaction cross-section 

= q{2J + ^f{m)M{s, t) (3.7) 

-\- background terms 

Averaging over t over an interval < ^1,^2 > gives a shape formula for the mass 
distribution 

/(.,m) = q{2J + l){-f^^^^^^—^f{m)N{s) (3.8) 

q (m^ — mj^y + mj^F^ 

+ background terms 
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where 

^(^) = / •^(^'^)^^= / \Qi'^t)\'dt (3.9) 

tl tl 

Setting /(m) = 1 and ignoring the background we get the Pisut-Roos resonance 
shape formula^^ for t-averaged mass distribution 

/(m) = NqF{m)\BW\'^ (3.10) 

where N is the normalization constant, q is the phase space factor, F{m) is the 
Pisut-Roos shape factor 



F(m)H2J+l)A^ = i^^ (3.11) 



and BW is the Breit-Wigner amplitude 



BW = , "^f . ^ (3.12) 

The Pisut-Roos resonance shape formula (3.10) has been extensively used to fit par- 
tial wave intensities in previous amplitude analyses of ttN — > tt'^tt~N on polarized 
targets (Ref. 14, 15, 16 and 23). 

B. Phenomenological resonance shape formula. 

In general, the experimental distribution I{m) in a certain mass region is fitted 
to a functional form^*^ 

I{m) = aRlR{m,mR,r) + aBlB{m) (3.13) 

where an and as give the fractions of resonant contribution and incoherent back- 
ground. Normally Ir is taken as a square of the Breit-Wigner amplitude multiplied 
by a phase space factor. A coherent term may be added to the Breit-Wigner ampli- 
tude, typically a constant term. In general, the background /^(/n) is a polynomial. 

In the case of n~p n~7T~^n reaction, the relevant phase space factor is just the 
pion momentum q in the tt+tt" cm. system and one can write for mass distributions 
in this reaction a phenomenological resonance shape formula 

I{m) ^ Nq{m){\BW\'^ + B} (3.14) 
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where N is overall normalization factor and B is the background term. We can 
take B = or B =constant. When 5 = 0, the phenomenological shape formula 
(3.14) is obtained from Pisut-Roos resonance shape formula (3.10) by setting their 
shape factor F = 1. We see from (3.11) that Pisut-Roos formula (3.10) converges 
to phenomenological formula (3.14) for large m when background B = 0. 

IV. The mass and width of cr(750) from fits to 5'-wave amplitude jSpE. 

As seen in Fig. 1 and 2, the Monte Carlo method and the method yield 
very similar results for the 5'-wave amplitudes |5'pE and |5'pE in n~p n~n^n 
at 17.2 GeV/c and for -t = 0.005 - 0.20 (GeV/c)^. Both methods show that 
the amplitude |S'pE resonates in both solutions while the amplitude |5'pE is non- 
resonating in both solutions. The Monte Carlo results appear to be smoother than 
the results. The Monte Carlo method found no physical solution at mass bin 
890 MeV. The solution found by method at this mass is far off from the general 
trend of data in solution 1 for IS'pE. For these reasons the mass bin 880-900 MeV 
was excluded from the fits to j/SpE. 

To determine the best values of the mass and width of cr(750) state from the 
mass distribution of the resonating amplitude |5'pE we used 4 types of fitting ap- 
proaches and used a criterion to determine the best fits. In the first approach 
we used a single Breit-Wigner fit. In the second approach we added an incoherent 
constant background to the single Breit-Wigner. In the third and fourth approaches 
we used two different versions of constant coherent background. In each approach 
we used both Pisut-Roos and phenomenological resonance shape formula and found 
they give very similar results. The inclusion of background leads to the narrowing 
of the width of (t(750). The best solution is obtained by the fourth approach 
leading to a conclusion that a(750) is a narrow state with a width about 100 MeV. 
The fitting was done using the CERN optimization program FUMILI.'^^ 

A. Single Breit-Wigner fit. 

In this approach the mass distribution |S'pE is fitted to a formula 



SfE = qFNs\BW 



(4.1) 




(4.2a) 
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for Pisut-Roos shape formula or 

F = 1. (4.26) 
for the phenomenofogical shape formula. BW is the Breit-Wigner amplitude 

BW = , "^f (4.3) 

mj^ — — zm^r 

where mu is the resonant mass. The mass dependent width T{m) depends on spin 
J and has a general form 

r = r,(^)"«^f(i«I> (4.4) 

qr Dj{qr) 

In (4.4) qji = q{m = mji) and Dj are the centrifugal barrier functions of Blatt and 
Weishopf:^^ 

Doiqr) = 1.0 

Di{qr) = 1.0 + (qrf (4.5) 

where r is the interaction radius. 

The results of the fit are shown in Fig. 3 and 4 for the Pisut-Roos and phe- 
nomenological shape formulas, respectively. The corresponding curves for both 
shape formulas are nearly identical. The numerical results are presented in Table 1. 
The fits to \S\^T, obtained by method have significantly higher values of x^/d.o.f. 
However both methods give a a mass in the range 730-750 MeV and a width in 
the range 230-250 MeV. Only the solution 1 of the method gives a lower width 
around 190 MeV. 

An important feature of the fits to \S\'^'E with single Breit-Wigner formula 
noticeable in Fig. 3 and 4 is that all fits lie below the maximum values of the 
mass distributions for each solution and the method of analysis. This inability of 
the single Breit-Wigner formula to reproduce the resonant shape of the amplitude 
\S\'^Tj suggests that background contributions are important and their effect on the 
mass and width of a state should be investigated, at least approximatively. 

B. Breit-Wigner fit with incoherent background. 

In this case we fit the mass distribution for |5'pS to a formula 

\S\'^E^qFNs{\BWf + B} (4.6) 
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where B is the incoherent background added to the Breit-Wigner formula (4.1). In 
general, S is a polynomial in m. However, since we have only 14 data points in the 
resonant mass range of 600-880 MeV, we will take B =constant. 

The results of the fit are shown in Fig. 5 for the phenomenological shape formula 
{F = 1). The results with Pisut-Roos shape formula are very similar. The numerical 
results are given in Table 2. We notice a dramatic improvement of the fit to the 
solution 2 for both methods which yields a better x^/dof and a narrow width of 
about 100 MeV. There is also some improvement of the fit to the solution 1 in 
particular for the method solution. This improvement is again associated with a 
lower x^/dof and a narrower width of 202 MeV and 147 MeV for the Monte Carlo 
and methods, respectively. The mass of the a state remains in the range of 
730-745 MeV. 

While the fits to solutions 2 are now much improved, the fits to solutions 1 
are still not satisfactory with the fitted curves still below the maximum values of 
these mass distributions. To make further progress we turn to coherent background 
contributions. 

C. Breit-Wigner fit with coherent background. 

The nonresonant behaviour of the amplitude |S'pS (recoil nucleon transversity 
down) strongly suggest the presence of a coherent nonresonating background. A 
part of coherent background also comes from the contribution of isospin 1 = 2 
amplitudes (see eq. (3.6)) which we neglected in the single Breit-Wigner fit. To 
understand the origins of the coherent background and to discuss its form for fits 
to liSpE it is useful to express the unnormalized moduli of S'-wave transversity 
amplitudes in terms of unnormalized helicity amplitudes. Using (2.13) we write 

\S\^E = l\So + iS\''E = qF\Fo + iF^\'^ 

|5|2e= ^|5o-i<Si|2E = gF|Fo-iFip (4.7) 

where -Fq and F\ are unnormalized S'-wave helicity amplitudes. The terms qF have 
the same meaning as in (4.1) and anticipate the use of (4.7) for Breit-Wigner fits 
to mass distribution of jiSpE. Near the resonance with mass we assume the 
following form of the helicity amplitudes 

F„(s, t, m) = i?„(s, t, m)BW{m) + S„(s, t, m) (4.8) 
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where n = 0, 1 is the nucleon hehcity flip, BW is the Breit-Wigner amphtude (4.3), 
Rnis, t, m) is the pole term and Bn{s, t, m) is the nonresonating background which 
includes the contribution from the nonresonating isospin 1 = 2 amplitudes. The 
energy variable s is fixed and will be omitted in the following. Since the experimental 
mass distributions are averaged over broad t-bins, we will eventually average also 
over the momentum transfer variable t. With the notation e = ±1, we can then 
write (4.7) in a compact form as follows 

Fo + ieFi = R^{t, m)BW{m) + B^{t, m) (4.9) 

where Re = Ro + leR, and Be = Bo + leB. It is useful to factor out the phase of R^ 
and define 

Re = \Re\e"^' 

Ce = Bee-"t>^ (4.10) 

Then (4.9) takes the form 

Fo + ieFi = {\Re\BW + CJe*"^^ (4.11) 
and the moduli squared of (4.7) read 

|Fo + ieFi|2 = + (ReCe)^ + (ImCe)^+ 

+2\Re\{ReCeReBW + linCelinBW} (4.12) 

We now recall that 

ReBW = ( ""^~"'' )|Wp = w\BW\'' 

ImBW=\BW\'^ (4.13) 

Hence 

|Fo + zeFi|2 = {\R,f + 2\R,\ReCeW + 2\R,\lmCe}\BW\'^ 

+(ReC,)2 + (ImC,)2 (4.14) 

Since the amplitude |5'pE(e = +1) does not show a clear resonant behaviour (Fig. 1 
and 2), we can conclude from (4.14) that the sum of terms in the parentheses must 
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be small or zero. This most likely means that \R+\ is small or zero implying that 
the pole terms in helicity amplitudes are related approximately as Rq ~ —iRi. 

For the resonating amplitude |5'pE(e = —1) the second and third terms in 
the parentheses in the equation (4.14) represent the effect of coherent background. 
In general the functions \R-\ and C_ will depend on both t and m. Since these 
functions are not known and since we have only 14 data points in the resonance mass 
region 600-880 MeV, we will work in the approximation of constant background. 
At this point there are two possibilities. 

(1) We assume that \R-\ and C- are constants independent of t and m. In this 
case no averaging over t is necessary and we can write (4.14) in the form 

\S\'^T. = qFNs{[l + 2wBi + 2B2] \BW\'^ + + S|} (4.15) 



where 



This possibility is equivalent to assuming that the constant parts of \R-\ and C_ 
dominate in the resonant mass range 600-880 MeV. We also notice that in this 
case the incoherent part B\ + i?| is correlated with the coherent contribution on 
(2wBi^-2B2)\BW\'^ in the formula (4.16) through the common parameters Bi and 



B 



2- 



(2) In the second possibility, we assume that \R-\ and C_ are both dependent on t 
and m. In this case we must average (4.14) over t over the experimentally measured 
interval < ti,t2 >• The averaging of (4.14) over t yields 

I^I^S = qF{[r + 2wa + 2h]\BW\^ + c} (4.17) 

where 

r =< >, a =< |i?_|ReC_ > 

b =< |i?_|ImC_ >, c =< (ReC_)2 + (ImC-)^ > (4.18) 

In (4.18) the symbol < > represents averaging over t over interval < ti,t2 >■ In 
general, the functions r, a, 6, c will depend on the mass m. Since we do not know 
these functions, we will assume constant values. But then there is no distinction 
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between r and 2b which can be combined into one parameter Ns = r + 2b as they 
are two constants in a sum. Then (4.17) has the form 



where Bi = b/Ns and B = c/Ns are the coherent and incoherent contributions 
to the resonance shape formula. This approximation is equivalent to assumption 
that the functions |-R_| and C_ depend mostly on t and only weakly on m. Notice 
that in this case the incoherent contribution B is not correlated with the coherent 
contribution as the parameters B and Bi are independent. 

We will refer to the first possibility (1) as Breit-Wigner fit with constant co- 
herent background and to the second possibility (2) as the Breit-Wigner fit with 
t-averaged constant coherent background. 

The results of the Breit-Wigner fit with constant coherent background are 
shown in Fig. 6 and Table 3. The results of the Breit-Wigner fit with the t-averaged 
constant coherent background are given in Fig. 7 and Table 4. Both Figures and 
Tables refer to the phenomenological shape formula with F = 1. The results with 
Pisut-Roos resonance shape formula (F given by (4.2a)) are very similar for the 
masses and widths although there are some differences in the fitted values of the 
constants Bi, B2 or Bi and B. 

An inspection of Figures 6 and 7 show much improved fits to the data on 
mass distribution of |S'pS. The overall best fit (as judged by the lowest values of 
X^/dof) is provided by the Breit-Wigner fit with the t-averaged constant coherent 
background. However the improvements in x^/dof appear only in solution 1 of 
Monte Carlo method and solution 2 of the method. Again, the Monte Carlo 
method achieves beter values of x^/dof compared to the x^ method of amplitude 
analysis. 

The improvements in the fits brought about by the inclusion of coherent back- 
ground have important consequences for the fitted values of the mass and width of 
cr(750) state. From Tables 3 and 4 we find that the mass of a in solution 1 is about 
30 MeV higher than the a mass found in solution 2. The Monte Carlo method gives 
the best value of a mass 774 MeV in solution 1 and 744 MeV in solution 2 (Table 
4). The x^ method gives the best value of a mass 761 MeV in solution 1 and 733 
MeV in solution 2 (Table 4). The data on polarized target cannot distinguish these 
two solutions. Since the two masses are close, we can work with a solution average. 




(4.19) 
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The solution average for a mass is 759 ± 22 MeV for Monte Carfo method and 747 
± 16 MeV for the method. The average over the two methods gives a mass 753 
± 19 MeV. 

The most significant effect of the inclusion of coherent background is the re- 
duction of the value of the width of a. The Monte Carlo method gives for the 
best value of a width similar values of 101 MeV and 103 MeV in solution 1 and 2, 
respectively (Table 4). The method gives for the best value of a width 134 MeV 
in solution 1 and 93 McV in solution 2 (Table 4). The data on polarized target 
cannot distinguish these two solutions, but the high values of x^/dof for method 
tend to favour the values for a width from the Monte Carlo method which has low 
values of x^/dof. The solution average for the a width is 102 ± 61 MeV for Monte 
Carlo method. The solution average for the a width is 113 ± 44 MeV for the x^ 
method. Since the error on the a width is larger for the Monte Carlo method, the 
two results are essentially compatible. The average over the two methods gives a 
width 108 ± 53 MeV. 

In conclusion, we propose to adopt the solution and method averages from the 
best fit values of Table 4 as the standard values of mass and width of the a state. 
The obtained values are 

m^ = 753±19MeV , = 108 ± 53 MeV (4.20) 

V. The interference with /o(980) in fits to amplitude |5'pE 

Tornqvist suggested^^ that the interference of (t(750) with /o(980) resonance 
could influence the determination of resonance parameters of a (750). In the old 
phase shift analyses (obtained using the invalid assumption of absence of Ai- 
exchange), the resonance /o(980) plays an important role of smoothly interpolating 
the "Down" solution for 5g below 900 MeV with the results for 5^ above 1000 MeV. 

We will now investigate the effect of interference of (7(750) with /o(980) on the 
determination of resonance parameters of cr(750). We will find that the effect is very 
small. This is consistent with the fact that /o(980) is a very narrow resonance and 
it is positioned sufficiently far away from the narrow and strong resonance cr(750). 

The experimental data in the /o(980) mass region are given in the CERN- 
Munich analysis-*^^ of tt~p Tr^Tr+n on polarized target at 17.2 GeV/c for dipion 
masses 600-1800 MeV. From Fig. 2 and Fig. 6 of Ref. 14 it is possible to reconstruct 
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the amplitudes \S\'^'E and The two solutions are shown in Fig. 8. The 

amplitude jSpS resonates at 750 MeV in solution 1 and at 800 MeV in solution 2. 
It shows a high value at 960 MeV and a pronounced dip at 1000 MeV, indicating 
an interference of /o(980) with background in this mass region around 1000 MeV. 
The structures are less dramatic in \S\'^T, which does not show a (750) but a dip at 
1000 MeV is still observable. 

To proceed, we extend our parametrization (4.8) of |S'pE to include /o(980) 
resonance. Recall from (4.7) that |<S'pS = qF\Fo — iFi\. Now we write for the 
helicity amplitudes Fq and Fi 

Fn = Fi''\s,t,m)BWa{m) + Rlf\s,t,m)BWf{m) + Bn{s,t,m) (5.1) 

where index a refers to (7(750) and / refers to /o(980). 
Then 

Fo - iFi = R^{s, t, m)BWa + i?/(s, t, m)BWf + B{s, t, m) (5.2) 

Assuming that the coefficients Ra,Rf and the background B are independent of t 
and m, we get an extension of the parametrization (4.15) 

{SfE = qFNs{[l + 2w^Bi + 2B2\\BW^\''+ (5.3) 

+Bf + Bl + [Cl + Cl]\BWf\^+ 
+2[{w^\BW^\'' + Bi){wfCi - C2)+ 
+{\BW^f + B2)iCi+WfC2)]\BWff} 

where 

WH='^^^^ , r = Tn{^) , R = aJ (5.4) 

If we assume that R^, Rf and B depend on t and perform t-averaging, the extension 
of parametrization (4.19) then reads 

\S\^J: = qFNs{[l + 2w^Bi] \BW^\^ + B+ 

+2[w^Bi + w^iwfCi - C2)\BWf\'^+ 
+{Ci + WfC2)\BWf\'']\BW^\''+ 

+ {Di+WfD2)\BWf\''} (5.5) 
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In the above parametrizations (5.3) and (5.5) the coefficients Ns, Bi, B2, (B), Ci, 
C2, -Di, D2 are real constants. The data between 900 and 1120 MeV exist only in 
40 MeV mass bins. Thus there is not enough data to fit the resonance parameters 
of /o(980). Instead we fix the mass of /o(980) at 980 MeV and its width at 48 MeV 
in the Breit-Wigner amplitude BWf. Also, in our fits we took for below 880 

MeV the results from our Monte Carlo analysis (in 20 MeV bins) and between 900 
and 1120 MeV we took the results of CERN- Munich analysis (in 40 MeV bins) from 
Fig. 8. 

The two parametrizations (5.3) and (5.5) yield virtually identical fits from 
600 to 1120 MeV and the same values for mass and width of a"(750). The fit for 
parametrization (5.3) is shown in Fig. 9 and the numerical values of the parameters 
are given in Table 5. There is a small improvement of x^/dof in Solution 1 which 
shows a better fit with the /o(980) interference. Comparison with the corresponding 
Table 3 shows a small increase in the mass of a in both solutions. There is a decrease 
of the a width in solution 1 from 114 MeV to 95 MeV and an increase in a width 
in Solution 2 from 104 MeV to 135 MeV. The solution averages are 

rria = 768 ± 22 MeV , = 115 ± 38 MeV (5.6) 

The effect of /o(980) interference is thus a small increase of average mass and width 
of a as compared to values in (4.20). It is not possible to claim'^^ that the low mass 
and the narrow width of a (750) are artifacts due to the neglect of interference of 
(7(750) with /o(980) in our fits. 

Both fits reproduce well the a resonance peaks below 880 McV in both solutions 
and the interference patterns between 920 and 1120 MeV. Particularly noteworthy 
in Fig. 9 is the dramatic drop in |5'pS between 960 and 1000 MeV due to destructive 
interference of /o(980) with the background. The good fit in this region suggests 
that the assumption of constant coherent background and resonance couplings is a 
good approximation. 

We have also attempted to fit the whole mass region of 600-1520 MeV us- 
ing a three resonance parametrization with a constant background and resonance 
couplings. The fit was not successful as the /o(1300) resonance was not well re- 
produced. This indicates that the background above 1120 MeV is different and the 
assumption of constant background for such a large mass range does not work. Next 
we fitted the /o(1300) resonance in the mass range of 1120 to 1520 MeV to a single 
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Breit-Wigner with incoherent background. The results are shown in Fig. 9 for the 
two solutions (differing in values of at 1480 MeV). The Solution 1 above 1120 

MeV connects smoothly with both solutions below 1120 MeV while the Solution 2 
shows a small discontinuity at 1120 MeV. Surprisingly, the incoherent background 
in both solutions is consistent with zero. This again indicates that above 1100-1200 
MeV the background (if any) is different from the low mass region below 1100 MeV. 
The numerical results of the fit to /o(1300) in the mass region 1120-1520 MeV are 
given in Table 6. We note the similarity of mass and width of resonances /o(1300) 
and /2(1270). 

VI. The mass and width of cr(750) state from the fits to ^-wave 
intensity Is- 

Previous amplitude analyses"^^"-^^'^^ of 7r~p — > Tr^Tr+n and Tr+n — > tt+tt^p 
data on polarized targets fitted only certain partial wave intensities using Pisut- 
Roos resonance shape formula without any background. It is of interest to perform 
Breit-Wigner fits to the S'-wave intensity Is and compare the results with the results 
of fits to resonating amplitude \S\'^'E in 7r~p 7r~7r+n at 17.2 GeV/c. Because of 
lower statistics, data for n^n 7r~^7:~p at 5.98 and 11.85 GeV/c allow fits only 
to S'-wave intensity Is- This is thus our primary aim in fitting 5"- wave intensity: 
to extract information about the mass and width of a in 7r~^n — > tt~^tt~p reaction 
measured at larger momentum transfers —t = 0.2 — 0.4 (GeV/c)^. 

Let us recall that the S'-wave intensity Is is defined as 

Is{s,t,m) = {\Sf + |5|2)S = {\Sof + (6.1) 

Since there are two independent solutions for the amplitudes \S\'^ and \S\'^, 
there are 4 solutions for the S'-wave intensity. We label these 4 solutions as 
Is{l, 1), /s(l, 2), /s(2, 1) and Is{2, 2) where 

Is{t,j) = i\S{t)f + \S{mE , ^,J = 1,2 (6.2) 

The results for the 4 solutions of Is obtained by the Monte Carlo amplitude analysis 
are shown in Fig. 10. The results for /^(l, 1) and /s(2, 2) obtained by the 
minimization method are shown in Fig. 1 1 . Again, there is a remarkable agreement 
between the results of these two different methods of analysis. The solutions /^(l, 1) 
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and 15(2, 1) are clearly resonating but the solutions /^(l, 2) and 15(2, 2) do not show 
a clear resonant behaviour. This is due to the large nonresonating contribution from 
the amplitude (see Fig. 1 and 2). The amplitude represents a nontrivial 

nonresonating background in all four solutions and is thus expected to distort the 
results of Breit-Wigner fits to Is- 

We first performed fits to Is using a single Breit-Wigner formula without any 
background 

Is ^qFNs\BW\'^ (6.3) 

In all fits to S'-wave intensities we used the Pisut-Roos shape factor F = (2J + 
l)(m/g)^. The results are shown as solid lines in Fig. 10 and 11 and in Tables 7 
and 8 for Monte Carlo and methods, respectively. We notice in Fig. 10 and 

11 that the single Breit-Wigner fit is well below the maximum values of the mass 
distribution Is- In Monte Carlo analysis the mass of a is around 766 MeV in all 
four solutions. The width is around 260 MeV for the first 3 solutions and is larger at 
303 MeV for the solutions 15(2, 2). In method the a mass and width in solution 
/^(l, 1) is in agreement with Monte Carlo results, but the width of /s(2, 2) is larger 
at 408 MeV and also mass is higher at 786 MeV. 

Next we performed Breit-Wigner fit with a constant incoherent background 
using a formula 

Is = qFNs{\BW\'^ + B} (6.4) 

where B is the constant background term. The results are shown as dashed lines 
in Fig. 10 and 11 and in Tables 7 and 8 for the Monte Carlo and methods, 
respectively. While the masses of a remain the same, there is a general reduction of 
the width of a associated with improved fits to the data and lower values of x^/dof. 
In Monte Carlo method the width of cr is reduced to 210 MeV in the first 3 solutions 
to Is- However the most dramatic and unexpected change occurs in the solution 
15(2,2) in both methods. There is a considerable improvement in the fit to the 
data and the width is drastically reduced to 188 MeV in both methods indicating 
the existence of a narrow a state even in the broad looking mass distribution. 

The best determination of a width from the fits to S-wave intensity Is is still 
double of the best value obtained in fits directly to the amplitude \S\'^T, (Table 
4). This discrepancy shows that the determination of resonance parameters from 
the spin-averaged intensities is not fully reliable when there is a presence of large 
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nonresonating nontrivial background as is the case of the amphtude j/SpS. The 
characteristic feature of this situation is that the S'-wave intensity does not show a 
clear reonant structure in all four solutions. 

This situation does not occur in the data on S-wsive intensity in Tr'^n — > 7r"'"7r~p 
at larger momentum transfers —t = 0.2 — 0.4 (GeV/c)^. The results from Monte 
Carlo amplitude analysis are shown in Fig. 12 and 13 at 5.98 and 11.85 GeV/c, 
respectively. We notice that all four solutions at both energies show clear resonant 
structures. This suggests that the determination of resonance parameters from S- 
wave intensities at these momentum transfers should be more reliable. However, 
this advantage is somewhat offset by the lower statistics of the data and large errors. 

We have again performed fits using single Breit-Wigner formula (6.3) and the 
Breit-Wigner fit with constant incoherent background using formula (6.4). The re- 
sults are shown in Fig. 12 and 13 and in Tables 9 and 10 for incident momenta of 5.98 
and 11.85 GeV/c, respectively. The fit with constant background (dashed lines) is a 
clear improvement over a single Breit-Wigner fit (solid lines). The improvement of 
the fit with the constant background is again associated with lower values of x^/dof 
and with reduction of the width of a in all solutions at both energies. However, 
there are differences in values for the mass and the width of a between the solutions 
as well as between energies. At 5.98 GeV/c, the mass ranges from 706 to 745 MeV 
and the width ranges from 145 to 262 MeV. At 11.85 GeV/c, the mass is higher 
and ranges from 756 to 782 MeV while the width is lower ranging from 117 to 202 
MeV. The differences are probably due to lower statistics. 

The solution averages for the mass and width of a from fits to Is are as follows: 
At 5.98 GeV/c 

rua = 730 MeV ± 27 MeV , = 195 ± 81 MeV (6.5) 
At 11.85 GeV/c 

= 768 ± 17 MeV , = 166 ± 54 MeV (6.6) 

At 17.2 GeV/c 

= 767 ± 9 MeV , F^ = 204 ± 75 MeV (6.7) 

The best values for the mass and width of a obtained from fits to the S-waye 
intensities at the three energies are in general agreement. The small differences are 
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likely due to the fact that the approximation of constant incoherent background 
may work differently at various energies and momentum transfers. The differences 
in mass of a from the fits to jS'pS and to Is are small. The difference in the value 
of the width from the best fits to |5'pS with coherent background and the fits to Is 
are somewhat large but the results are still consistent. At 17.2 GeV/c they are due 
to large nonresonating contributions from the amplitude The differences also 

reflect the need for inclusion of coherent background and its better description than 
a constant. This in turn would require more data of high statistics in the resonance 
region 600-900 MeV. 

VII. Remarks on determinations of tttt phase shifts. 

The amplitude analyses of measurements of ttNi; — > 'k'^'k~N on polarized tar- 
gets provide a model-independent and solution-independent evidence for a narrow 
scalar state 1 = 0"''"'"(750). The question arises how to understand the absence of 
such a state in the conventional /S-wave phase shift 5q in tttt scattering.^^'^^"^-*^ 

Of course, there are no actual measurements of pion-pion scattering and there 
is no partial-wave analysis of tttt — > tttt reactions in the usual sense. The tttt phase 
shifts are determined indirectly from measurements of tt~p TT~TT'^n on unpo- 
larized targets using several strong enabling assumptions. One of these crucial 
assumptions - the absence of Ai exchange amplitudes - leads to predictions for po- 
larized spin density matrix (SDM) elements and for the measured amplitudes, and 
it is thus directly testable in the measurements on polarized targets. As we shall see 
below, the assumption of absence of Ai-exchange amplitudes is totally invalidated 
by the data on polarized targets. The polarization measurements also cast some 
doubt on the fundamental assumption of factorization of mass m and momentum 
transfer t in the crucial pion exchange amplitudes. We must use the results of mea- 
surements on polarized targets to judge the validity of tttt phase shifts, and not vice 
versa. We are thus led to the conclusion that the indirect and model-dependent 
determinations of tttt phase shifts cannot be correct. This explains the absence of 
1 = 0"''"'"(750) resonance in the 5q phase shift from these analyses. 

We will now review the basic assumptions common to all determinations of tttt 
phase shifts. "'24-31 

A priori, there is no connection between the partial wave amplitudes in 
TTTT — > TTTT Scattering and the production amplitudes in t:N — > n~^n~N reactions. 
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We recall that in ttN — > tt+tt^A?" there are two 5'-wave production amplitudes 
S{s,m,t) and S{s,'m,t) (or So{s,m,t) and Si{s,m,t)) while in tttt — > tttt there is 
one S'-wave amplitude (or phase shift 5q) dependent only on the energy E. Also, 
in ttN — > 7r"'"7r~A'" there are six P-wave production amplitudes L, L, U, U, N, N (or 
Ln,Un, Nni n — 0, 1) which depend on variables s,m,t while in tttt — > tttt there is 
again one P-wave amplitude (or phase shift dl) dependent only on the energy E. 
To make the connection between the production amplitudes in ttN — > 7r"'"7r~A^ and 
the partial-wave amplitudes in tttt tttt the following assumptions of factoriza- 
tion and identification are postulated in all determinations of tttt phase shifts from 
unpolarized data on ttN — > tt~^tt~N. 

The starting point are the dimeson helicity A = pion exchange amplitudes 
Si and Li in the t-channel. It is assumed that the t and m dependence in these 
amplitudes f actor izes: 



-t ^ , , m 



Siism,t) = N^^Fo{t) — foim) 



Li{s,m,t) = N^^Fi{t) — h{m) (7.1) 

where t is the momentum transfer at the nucleon vertex, m and q are the dipion mass 
and the tt~ momentum in the tt~tt'^ cm. frame. The form factors fj{t) describe 
the t-dependence and the functions fj{m), J = 0, 1, describe the mass dependence. 
N is a, normalization constant. Furthermore, the functions fj{m) are assumed to 
be the partial-wave amplitudes in tt~tt'^ — > tt~tt'^ reaction at cm. energy m: 

Jo — g/o "T g/o 

fi = ft' (7.2) 

The partial wave amplitudes fj with definite isospin / are defined so that in the 

TTTT elastic region 

fj = sin S'je^''^ (7.3) 

The phase shifts 5j are determined from the amplitudes and Li which are calcu- 
lated from the data on tt~p —>■ tt^tt+tt, on unpolarized target. However the calcula- 
tion of amplitudes Si and Li from the tt~p TT~TT'^n data cannot be done without 
additional assumptions. There is simply more amplitudes than data. To proceed 
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further all determinations of tttt phase shifts must assume that all ^i-exchange 
amplitudes vanish: 

So = Lo = Uo = (7.4) 

With the assumptions (7.4), two solutions for the ^S-wave phase shift 6q are 
found:^^'^^ "Down" solution which is non-resonating and "Up" solution which res- 
onates at the mass around 770 MeV with a width about 150 MeV. The resonating 
solution was rejected because it disagreed with the 7r°7r° mass spectrum from a 
low-statistics experiments^ on 7r~p — > 7r°7r°n at 8 GeV/c. 

There is no theoretical proof of factorization (7.1) and identification (7.2) of 
functions fj with tttt partial-wave amplitudes. It is not obvious that the tttt phase 
shifts calculated from 7r~p — > 7r~7r+n data using the assumptions (7.1)-(7.3) would 
coincide with tttt phase shifts determined directly from real pion-pion scattering. 
Only such comparison could test the assumption (7.2). 

The factorization (7.1) implies that the mass spectrum of amplitudes \Si\'^ 
and \Li\'^ is independent of t. This consequence of factorization can be tested 
in measurements of nN 'k^'k~N on polarized targets. In Fig. 14 we show t- 
evolution of mass dependence of lower and upper bounds^*^ on normalized moduli 
|L|2, |L|2, and |t7p. The data at t = -0.068 (GeV/c)^ are from Ti'p ^ Tf-Tr+n 
at 17.2 GeV/c, the rest is from 7r+n n~^n~p at 5.98 GeV/c. The Fig. 14 shows 
a clear and pronounced dependence of mass spectra of amplitudes |Lp and |Lp 
on momentum transfer t. In particular, there is a clear change of mass spectrum 
below —t = 0.25 (GeV/c)^, a region of t relevant to determinations of tttt phase 
shifts. While this change could be entirely due to exchange amplitude Lq, this 
cannot be guaranteed. The factorization assumption (7.1) thus cannot be taken for 
granted and further tests of this assumption are required in future high statistics 
measurements of ttN tt+tt^AT on polarized targets. 

The assumption (7.4) of absence of Ai-exchange amplitudes has several conse- 
quences that can be directly tested in measurements on polarized targets. ^From 
(2.13) we see that absence of Ai exchange amplitudes implies 

|A| = |A| for A^S,L,U (7.5) 

The equality of moduli of amplitudes with the recoil nucleon transversity "down" 
and "up" is not observed experimentally. We can see in Fig. 1 and 2 that the .S-wave 
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amplitudes and are clearly unequal at 17.2 GeV/c and —t = 0.005 — 0.20 
(GeV/c)^. In Fig. 14 we see that the P-wave amplitudes \L\ and \L\ are different in 
every t-bin from 0.005 to 0.60 (GeV/c)^, and that the difference is largest at small 
t, the region of most importance to determination of tttt phase shifts. 

The Ai-exchange is large and nontrivial also above 900 MeV and in higher 
partial waves D and F. This finding of CERN-Munich analysis^^ of 
data on polarized target in the mass range 600-1800 MeV is shown in Fig. 15. 
The figure shows the ratios of moduli of amplitudes with recoil nucleon transversity 
"down" and "up" for S-, P-, D- and F-wave amplitudes with dimeson helicity A = 
which are directly relevant for determination of the corresponding phase shifts. The 
deviations from 1 indicate the strength of Ai-exchange. We can see in Fig. 15 that 
Ai-exchange is important in all waves up to 1800 MeV at small —t — 0.005 — 0.20 
(GeV/c)^. The determinations of tttt phase shifts above 900 MeV also assumed the 
absence of ^i-exchange amplitudes. We must conclude that the determinations of 
TTTT phase shifts from 5'- wave to F-wave in the mass region from 600 to 1800 MeV 
are not reliable. Theoretical calculations and analyses based on these phase shifts 
are therefore not reliable as well. 

Below 1000 MeV, where the S- and P-wave dominate, the assumptions (7.4) 
lead to predictions for polarized SDM elements that can be directly compared with 
the data. The predictions of (7.4) are^^ 

Pi + Plo + 2p^i = -2(p^o - Pli) = (7.6) 

Rep^o = Repf, = Rep^, = (7.7) 

The data for polarized SDM elements clearly rule out these predictions as is shown 
in Figs. 16 and 17 for ir'p tt-tt+ti at 17.2 GeV/c. We find that pf^ + p^g + 2pl^ 
and — 2(poo — Pii) have large magnitudes but opposite signs while 2p^_]^ has a small 
magnitude. The interference terms Rep^g, Repi^ and Repg^ are all dissimilar and 
have large nonzero values. On the basis of this evidence we again must conclude that 
the past determinations of tttt phase shifts from unpolarized data on ttN — > Tr+Tr^A?" 
are questionable. 

The assumption of absence of exchange amplitudes means that pion pro- 
duction in ttN — > TT+TT^AT reactions does not depend on nucleon spin. What the 
measurements of ttN Tr+Tr^iV on polarized targets found is that the pion produc- 
tion depends strongly on nucleon spin. The dynamics of the pion production is not 
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as simple as has been assumed in the past determinations of tttt phase shifts. New 
determinations of tttt phase shifts are now required that do take into account the 
existence of Ai exchange. Since the contributions of Ai exchange amplitudes are 
large and nontrivial, the revisions of tttt phase shifts will be significant. The new 
revised ^S-wave phase shift 5q is then expected to show evidence for narrow scalar 
state cr(750) in agreement with the measurements on polarized targets. 

VIII. Questions concerning evidence for narrow a (750) 
A. Up-Down ambiguity and analyticity constraints. 

Recently it has been claimed^^'^° that tttt phase shift 5q can be determined from 
the S-wave intensities Is obtained in our amplitude analysis of 7r~p Tr^Tr+n on 
polarized target at 17.2 GeV/c, and that it would show the old Up-Down ambiguity 
of 6q. Only Up solution indicates a narrow a state and it is excluded because it is 
inconsistent with Roy equations. From this it was concluded that a"(750) does not 
exist or that the evidence must be treated with reservation.^*^ 

To answer this objection we first recall from (2.15) that 

/s = (|^o|' + |^i|')S (7.1) 

Here the amplitude Si is connected to Sq through (7.1) and (7.3), and 5*0 is the 
unknown Ai exchange amplitude. It is obvious from this expression that the de- 
termination of ^0 from data on Is depends on the model used for Ai exchange 
amplitude Sq. The data on polarized target require large Ai exchange amplitudes. 
At present the Ai exchange amplitudes are not known. We must therefore conclude 
that the phase shift Sq cannot be determined from the data on S-waive intensity Is 
at present. 

Nevertheless, the data on Is do tell us something very important about the 
solutions for 5q. There are four solutions for /g: 15(1, 1 ),..., 15(2, 2). Consequently 
there will be a fourfold ambiguity in Sq for any given model of Ai exchange am- 
plitude ^o. However, as can be seen in Fig. 10, the four solutions for Is are all 
very similar quantitatively. Consequently the four solutions for Sq are expected to 
be very close to each other and similar. This contrasts with the large differences 
between the old Up and Down solutions. Fig. 18 shows S-wave intensity normal- 
ized to 1 at maximum for Down (curve A) and Up (curve B) solutions from the 
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typical analysis of Estabrooks et al. "'^ The large differences between the Up and 
Down solutions contrast sharply with the small differences shown between S-wave 
intensities 15(1, 1) and 15(2,2) in Fig. 18. On the basis of the similar behaviour 
of all solutions for Is we do not anticipate the emergence of the old Up-Down am- 
biguity problem in Sq. It is even possible that the small differences between the 
four solutions for Is can be explained entirely as a small ambiguity in Ai exchange 
amplitude leading to a unique determination of Sq from the data on polarized 
target. 

The above discussion applies also to the determination of P-wave phase shift 
dl from II — (|-^oP + The amplitude Li is connected to 6l by (7.1) and 

(7.3) while Lq is another unknown ^i-exchange amplitude. The four solutions for 
II are again very close so we expect similar solutions for 5}. 

Assuming a model for Ai exchange amplitudes and Lq, the obtained phase 
shifts 5q and 5l can be tested for consistency with dispersion relations^^ (Roy equa- 
tions). If an inconsistency is found it means that we have to modify our model for 
Ai exchange amplitudes Sq and Lq, and try again. It is important to realize that 
Roy equations do not test the validity of the experimentally measured amplitudes 
\S\^, \S\'^, |Lp, |Lp or intensities Is and II- The Roy equations are constraints 
only on tttt phase shifts which follow from the analyticity properties of partial wave 
amplitudes in tttt — > tttt scattering. However the requirement of consistency of phase 
shifts with the Roy equations can be used to constrain the possible models of Ai 
exchange amplitudes. 

We conclude that the experimental evidence for the narrow state a (750) is 
not in contradiction with analyticity and dispersion relations for nn partial waves. 
The existence of Ai-exchangc and narrow cr(750) are experimental findings from 
measurements on polarized targets independent of the Roy equations. These exper- 
imental facts cannot be refuted by comparisons with standard phase shifts because 
these were obtained using an invalid assumption of absence of ^i-exchange. 

B. Absence of (t(750) in 77 — > tt+tt" and central production pp — > ppTr+Tr". 

Morgan and Pennington suggested to discount the evidence for existence of 
narrow cr(750) in ttN — > tt+tt^A?" because this state has not been observed in 77 
7r"''7r~ reaction'*^ and in central production^^'^^ pp pp7r+7r~. However there are 
good reasons why one would not expect to observe narrow a(750) in these processes. 
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In the next section we shall argue that the narrow cr(750) is the lowest mass 
scalar gluonium 0"''"''(^^). The principal support for this proposal is precisely the 
fact that cr(750) state is not observed in 77 — > 7r"''7r~ reaction. Since gluons do 
not couple directly to the photons, we expect cr(750) not to appear in reaction 
77 — > 7r"''7r~ if it is pure gluonium or if it has only a small qq component. 

The reaction pp — pptt~^tt~ was measured^^ at the CERN Intersecting Storage 
Rings (ISR) in a search for scalar gluonium. The structures reported in the moments 
H{11) and H{31) near m(7r+7r") ^ 750 MeV are consistent with ct(750) and p°(770) 
interference. 

Assuming parity conservation there are 5 5'-wave amplitudes and 15 P-wave 
amplitudes in this reaction. The cr(750) state may contribute only to some 5'- wave 
amplitudes and not to the others, as it does in 7r~p 7r~7r+n with amplitudes 
and As we see in Fig. 10, the S'-wave intensity 15(2, 2) does not immediately 

suggest existence of narrow cr(750). With 5 S-waye amplitudes in pp — > pp7r'^7r~ it 
is very likely that cr(750) stays hidden. We can observe cr(750) in 7r~p — > tt^tt+tt- 
and TT'^n — > 7r"''7r~p reactions only when these production processes are measured on 
polarized targets, and the S- and P-wave amplitudes can be separated in a model 
independent way. For the same reasons we may see cr(750) in central production 
pp — > ppTT'^TT~ only when measurements with polarized initial protons are made and 
the resonating S'-wave amplitudes can be isolated. The ISR experiment does not 
separate the S- and P-wave amplitudes, and thus it is not conclusive. 

C. Comparison with other results for a state. 

DM2 Collaboration measured^^ 7r"'"7r~ mass distribution in J/^ip wtt+tt" 
decays and observed a quite broad low mass resonance (see Fig. 13a of Ref. 54). 
Interpreted as an / = O""""*" a state, a single Breit-Wigner fit gives = (414 ±20) 
MeV, Fo- = (494 ± 58) MeV. There is no indication for such state in our data on 
S-wave intensity Is in Tr+n — > tt+tt^p at 5.98 and 11.85 GeV/c (see Fig. 12 and 13 
above). The reasons for the discrepancy are not clear at the present. 

Several recent theoretical analyses^^'^^'^^ claimed existence of a cr meson with 
a mass around 1000 MeV and a broad width of 460-880 MeV. These analyses use as 
an input the /S-wave phase shift Sq and thus neglect the Ai exchange and other spin 
effects observed in pion production (see e.g. eq. (5) in Ref. 56). It is possible that 
when these analyses include in their fits Ai exchange that they will find a narrow 
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cr in agreement with the CERN data on polarized targets. 



IX. Constituent structure of the (t(750) resonance. 

In the usual quark model meson resonances are qq states. The mass of cr(750) 
is too low for it to be a qq state. The mass M of the qq state increases with its 
angular momentum L as M = Mo(2n+L) where n is the degree of radial excitation. 
The lowest mass scalar mesons are ^Pq states with masses expected to be around 
1000 MeV or higher. 

It was suggested that 0++(700) could be a four-quark qqqq state in the MIT 
bag model. However, more detailed studies of qqqq systems conclude that pure 
multiquark hadrons do not exist^^'^^ with tt+tt" decay. We can also exclude the 
possibility that (t(750) is a hybrid state qqg. The lowest mass hybrid state must 
be a or 1 state. Calculations based on bag models, QCD sum rules, lattice 
QCD and a string model all estimate^^ the masses of 0~^~^{qqg) states to be above 
1500 MeV. 

Ellis and Lanik discussed the couplings of scalar gluonium a on the basis of the 
low energy theorems of broken chiral symmetry and scale invariance, implemented 
using a phenomenological lagrangian.^^ They obtained for a — > 7r"'"7r~ decay the 
following partial width 

where Go =< 0\{as/Tr)F^i,F^'^\0 > is the gluon-condensate term^^ parametrizing 
the non-perturbative effects in QCD. The numerical values were estimated by the 
ITEP group64 be Go w 0.012 (GeV)^ or up to Go ^ 0.030 (GeV)^ in later 
calculations. ^^'^^ It is very interesting to note, that when we take Go = 0.015 
(GeV)^ the EUis-Lanik theorem (9.1) predicts partial width of cr ^ tt+tt" decay 
r = 107 MeV for the mass = 753 MeV. This result is in perfect agreement with 
(4.20), the solution and method average values of mass and width of cr(750) from 
the best fit to the measured mass distribution |5pS (Table 4). When we use for 
nia the value 768 MeV obtained in interference fits with /o(980) then EUis-Lanik 
theorem predicts a width T{a — > tt+tt") = 118 MeV, again in perfect agreement 
with (5.6) where r^. = 115 ± 38 MeV. From this agreement we can conclude that 
the cr(750) is best understood as the lowest mass gluonium state 0~^~^{gg). 

The gluonium interpretation of a (750) gathers further support from the lack 
of observation of cr(750) in the reactions 77 — > tt+tt" and 77 — > tt^tt^. Since gluons 
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do not couple directly to photons we expect cr(750) not to appear in reactions 
77 — > TTTT if it is a pure gluonium state or if it contains only a small qq component. 
This conclusion is supported by the PLUTO and DELCO data.^^'^® However, the 
more recent DMl/2 data^^'^° show an excess over the Born term expectation that 
is attributed to the formation of a broad scalar resonance with a two-photon width 
of (10 ± 6) MeV. This would suggest some qq component in the cr(750) state. The 
most recent results^^ are on 77 — tt^tt*^ which show no evidence for a scalar state 
near 750 MeV. 

Lattice QCD calculations by several groups^^"^^ initially concluded that the 
gluonium ground state 0~^~^{gg) has a mass near the meson: 740± 40 MeV. 
The most recent lattice QCD calculations predict much higher mass of the lowest 
scalar gluonium: the UKQCD group^^ predicts 1550 ± 50 MeV while the IBM 
group^^'^^ predicts 1740 ± 70 MeV. However, it is important to remember that 
these calculations are for quenched QCD so that there is no coupling of the primitive 
gluonium to quarks. The coupling of gluonium to two pseudoscalars may have a 
significant effect on the gluonium mass and width. 

We conclude that while the gluonium interpretation of the cr(750) state is in 
agreement with low energy theorems of broken chiral symmetry and scale invariance, 
it is at variance with the most recent lattice QCD calculations. It is necessary to 
study this discrepancy and understand its origins and implications. 

Finally we note that the anomalous energy dependence of pp and np elastic po- 
larizations and the departure from the mirror symmetry in ttN elastic polarizations 
at intermediate energies require a low- lying Regge trajectory^^''^'^ corresponding to 
a"(750). These anomalous structures in the polarization data may have been the 
first evidence for a gluonium exchange in two-body reactions. 

X. Amplitude spectroscopy — a new direction in hadron spectroscopy. 

The vast majority of hadron resonances have been identified through study of 
mass distributions of spin averaged cross-sections. Experiments with polarized tar- 
gets opened a whole new approach to experimental hadron spectroscopy by making 
accessible the study of hadron production on the level of spin dependent production 
amplitudes. We may refer to this new approach to detecting and studying hadron 
resonances as amplitude spectroscopy. 

This work represents the first effort to determine resonance parameters directly 
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from a measured spin dependent production amplitude. In this case it was the 
amphtude measured in 7r~pi — > Tr^Tr+n at 17.2 GeV/c. We have found that 

the best Breit-Wigner fits to the resonance cr(750) on the amphtude level differ 
markedly from the Breit-Wigner fits to the spin-averaged 5'-wave intensities, and 
provide a more reliable information about the resonance parameters. 

However, the significance of amplitude spectroscopy goes far beyond the more 
precise determinations of resonance parameters. The amplitude spectroscopy opens 
several prospects for new physics: 

(a) Let us call dominant resonances those states which can be experimentally iden- 
tified in spin-averaged cross-sections or Dalitz plots. Historically these dominant 
resonances led to the standard quark model and to QCD. However, a new species of 
subdominant resonances may exist which can be identified only at the level of spin 
dependent production amplitudes. The existence and properties of subdominant 
resonances could reveal new components of hadron structure to which unpolarized 
experiments are totally blind and could lead us beyond the standard quark model 
and standard QCD. 

(b) The production of resonances (dominant and subdominant) may depend on 
nucleon spin and this dependence will provide important information about the 
dynamics of hadron interactions and about the very nature of hadron resonances. 

(c) Standard QCD predicts new kinds of resonances such as dibaryon, gluonium 
and hybrid states. Many of these states may not be observable in the spin-averaged 
measurements which could explain the limitted success in identifying these states 
so far. 

The (t(750) is the first example of a subdominant resonance observable only on 
the level of measured spin dependent production amplitudes. We have interpreted 
this resonance as the lowest mass gluonium 0"''"'" (gg) in the Section IX. However, we 
must be open also to the possibility that cr(750) represents the first signal of a new 
physics beyond the standard quark model and QCD. It may indicate the existence 
of a new component in hadron structure. 

The measurements^® of K'^n — > K'^7r~p reactions at 5.98 and 11.85 GeV/c at 
CERN-PS also allowed a model independent amplitude analysis®-*^ of this reaction 
at 5.98 GeV/c. The results of the amplitude analysis are in excellent agreement 
with Additive Quark Model predictions.®^ The data also suggest^®'®^ the existence 
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of another subdominant resonance, I = \ 0"'""''(890) with a narrow width of about 
20 MeV. Such new resonance k(890) under i^°*(892) could also signal new physics 
beyond the standard QCD. One such possibility is discussed in Ref. 81. 

We recall that complete measurements of spin observables in two-body and 
quasi two-body reactions enable construction of spin amplitudes. The spin ampli- 
tudes show dip structures in their moduli associated with rapid and large changes 
of their relative phases. These dip structures resemble absorption resonances and 
systematic study of dip structures in two-body and other exclusive processes like 
ttN — > 7T~^7T~N is a natural extension of hadron spectroscopy into the space-like 
region. Study of time-like resonances and space-like dips in spin dependent ampli- 
tudes should bring entirely new insights into the hadron dynamics and structure. 

To explore these new frontiers of hadron spectroscopy and hadron dynamics, 
new advanced hadron facilities dedicated in large part to measurements with spin 
will be required. The proposed Canadian KAON Factory,^^ Los Alamos Hadron 
Facility®^ and European Hadron Facility*^ could integrate the new advanced tech- 
nologies of polarized beams and targets in a single spin physics facility^^ to system- 
atically advance the exploration and development of these new frontiers in hadron 
physics.^^"^^ 

XI. Summary 

The measurements of reactions 7r~p| Tc~n~^n at 17.2 GeV/c and yr+nj 
n~^TT~p at 5.98 and 11.85 GeV/c on polarized target provide model-independent 
and solution-independent evidence for a narrow scalar state a {750). The amplitude 
analyses of 7r~p^ — > 7r~7r^n at small t using minimization method^"^ and Monte 
Carlo method^^ yield very similar results for moduli of transversity amplitudes and 
cosines of their relative phases. In particular they agree that the transversity "up" 
S-wave amplitude resonantes near 750 MeV while the transversity "down" 

amplitude j^pS is nonresonating and constitutes a large background in the spin- 
averaged S-wave intensity Is = {\S\'^ + |5'p)E. For this reason it is preferable to 
determine resonance parameters of cr(750) directly from the measured mass distri- 
bution of |5pS. 

We have performed several types of Breit-Wigner fits to jS'pE. We have shown 
that the Pisut-Roos resonance shape formula and phenomenological shape formula 
give similar results. Single Breit-Wigner fits yield a width of cr (750) in the range 
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192-256 MeV. We have studied the effect of background in three approaches: inco- 
herent background, constant coherent background, and t-averaged constant coher- 
ent background. The last method yields the best fit with the lowest x^/dof. The 
solution and method average for the a mass and width from this best fit are 

m^ = 753±19MeV , = 108 ± 53 MeV (11.1) 

We also performed the conventional fits to spin-averaged 5- wave intensity Is- 
We found again that the inclusion of background (incoherent in this case) reduces 
the fitted value of the a width and improves x^/dof. Nevertheless, the direct fits to 
are preferable at 17.2 GeV. Due to lower statistics at 5.98 and 11.85 GeV/c, we 
must use results for Is to obtain a resonance parameters. All four solutions resonate 
at these larger momentum transfers but yield a broader a width: r^. = 195 ± 81 
MeV at 5.98 GeV/c and = 166 ± 54 MeV at 11.85 GeV/c. 

We conclude that the best overall estimate of the mass and width of cr(750) 
are the values in (11.1) from the best fit to \S\'^T, (Table 4). 

We have also examined the interference of (t(750) with /o(980) and found that 
it has only a small effect on the mass and width of (t(750). A fit to amplitude \S\'^T, 
in the mass range above 1120 MeV shows evidence for a scalar state with average 
mass 1280 ± 12 MeV and width 192 ± 26 MeV. 

The conventional 5'-wave phase shifts Sq show no evidence for the narrow 
cr( 750) state. It must be reiterated, that the past determinations of nn phase shifts 
from unpolarized data on TT~p —>■ n~Tv^n assumed the absence of Ai-exchange 
amplitudes. This assumption is invalidated by measurements of tt~p tt^tt+tt-, 
TT'^n — > 7T~^7T~p and K^n K~^n~p on polarized targets which find large and 
nontrivial ^i-exchange contributions. New determinations of tttt phase shifts are 
required that do take into account the existence of ^i-exchange. Since ^i-exchange 
contributions are large, the revisions of tttt phase shifts will be significant and should 
provide evidence for a narrow cr(750) state in agreement with the CERN data on 
polarized targets. 

The mass of cr(750) is too low for it to be a qq state. We proposed to identify 
cr(750) with the lowest mass scalar gluonium 0"''"''(^^). This proposal is supported 
by the perfect agreement with the EUis-Lanik theorem (9.1) relating the decay width 
of scalar gluonium T{a ^ 7r"''7r~) to its mass m^. Another experimental support 
for the gluonium interpretation of a(750) is its absence in 77 — > tt+tt" reaction. 
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However, the low mass of cr(750) is at variance with the more recent calculations of 
lattice QCD which predict masses of scalar gluonium above 1500 MeV. 

Experiments with polarized targets have opened a whole new approach to ex- 
perimental hadron spectroscopy by making accessible the study of hadron produc- 
tion on the level of production spin amplitudes. We may expect that this new field 
of amplitude spectroscopy will be further developed at the new proposed advanced 
hadron facilities. ^"^"^^ 
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xVdof 

Solution (MeV) (MeV) 



Pisut-Roos shape formula 



1(MC) 


736 ± 6 


230 ± 32 


1.40 ± 0.12 


0.388 


2(MC) 


745 ± 12 


240 ± 59 


1.71 ± 0.23 


0.276 




738 ± 4 


191 ± lb 


1.50 ± 0.11 


0.662 




752 ± 10 


253 ± 46 


1.79 ± 0.15 


0.968 


Phenomenological shape 


formula 






1(MC) 


732 ± 6 


231 ± 33 


6.50 ± 0.57 


0.418 


2(MC) 


740 ± 11 


241 ± 61 


7.94 ± 1.11 


0.288 




733 ± 4 


192 ± 16 


7.00 ± 0.50 


0.740 




747 ± 10 


256 ± 47 


8.29 ± 0.74 


0.986 



Table 1. Results of the fits to the mass distribution jS'pE measured in 7r~p — > 
TT^TT+n at 17.2 GeV/c using a single Breit-Wigner formula (4.1). The notation MC 
and indicates the solutions obtained by the Monte Carlo and minimization 
methods, respectively. 



48 



Solution (MeV) (MeV) 
Phenomenological shape formula 
1(MC) 731 ± 6 202 ± 110 
2(MC) 744±14 103 ± 74 
I(X^) 736 ±4 147 ± 43 

2(x^) 745 ±41 98 ± 41 



B Ns xVdof 



0.15 ± 0.55 5.75 ± 2.54 0.416 

0.73 ± 0.49 5.11 ± 1.84 0.144 

0.19 ± 0.17 6.13 ± 0.84 0.696 

0.70 ± 0.28 5.70 ± 1.30 0.626 



Table 2. Results of the fits to the mass distribution j/SpE measured in 7r~p — > 
TT^TT+n at 17.2 GeV/c using a Breit-Wigner formula with a constant incoherent 
background (4.6). The notation MC and as in Table 1. 
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Si B2 Ns xVdof 

Solution (MeV) (MeV) 
Phenomenological shape formula 



1(MC) 


770 ± 19 


114 ± 17 


0.90 ± 0.43 


1.34 ± 0.84 


1.09 ± 0.61 


0.136 


2(MC) 


745 ± 31 


104 ± 76 


0.02 ± 1.07 


1.84 ± 1.14 


1.09 ± 0.90 


0.144 


i(x^) 


761 ± 13 


138 ± 19 


0.34 ± 0.16 


0.69 ± .045 


2.41 ± 1.14 


0.362 




738 ± 20 


103 ± 112 


-0.17 ±0.65 


1.30 ± 0.50 


1.79 ± 0.86 


0.898 



Table 3. Results of the fits to the mass distribution jS'pE measured in 7r~p — > 
TT^TT+n at 17.2 GeV/c using a Breit-Wigner formula with constant coherent back- 
ground (4.15). The notation MC and as in Table 1. 
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rria B B Ng X^/dof 

Solution (MeV) (MeV) 

Phenomenological shape formula 



1(MC) 


774 ± 14 


101 ± 44 


0.73 ± 0.31 


0.29 ± 0.14 


3.99 ± 0.96 


0.108 


2(MC) 


744 ± 31 


103 ± 79 


0.73 ± 0.54 


0.01 ± 0.23 


5.10 ± 1.90 


0.144 


i(x^) 


761 ± 12 


134 ± 41 


0.25 ± 0.17 


0.15 ± 0.07 


5.74 ± 0.82 


0.362 




733 ± 20 


93 ± 48 


0.80 ± 0.39 


-0.12 ± 0.19 


5.31 ± 1.53 


0.592 



Table 4. Results of the fits to the mass distribution |S'pS measured in 7r~p — > 
TT^TT+n at 17.2 GeV/c using a Breit-Wigner formula with t-averaged constant co- 
herent background (4.17). The notation MC and ^ in Table 1. 
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l-SpS 5i B2 Ns xVdof 

Solution (MeV) (MeV) 

Phenomenological shape formula 

1(MC, x^) 778 ± 13 95 ± 27 1.20 ± 0.35 0.85 ± 0.33 1.24 ± 0.39 0.096 

2(MC, x^) 758 ± 32 135 ± 49 0.54 ± 0.69 1.00 ± 0.28 1.85 ± 0.83 0.162 

Ci C2 

1(MC, x^) 0.42 ± 0.52 1.25 ± 0.39 

2(MC, x^) -0.35 ± 0.67 0.97 ± 0.55 



Table 5. Results of the fit to the mass distribution l^pS in the mass range from 
600 to 1120 MeV taking into account the interference of cr(750) with /o(980) using 
the parametrization (5.3). The notation MC and x^ as in Table 1. 



l^pS m V B N xVdof 

Solution (MeV) (MeV) 

I(X^) 1284 ± 12 209 ± 29 0.001 ± 0.32 5.96 ± 0.62 1.393 

2(x^) 1276 ± 11 175 ± 24 0.001 ± 0.09 6.21 ± 0.70 1.738 



Table 6. The results of the fit to the mass distribution [S'pS in the /o(1300) mass 
region from 1120 to 1520 MeV using a single Breit-Wigner formula with incoherent 
constant background. The notation x^ as in Table 1. 
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Is B Ns xVdof 

Solution (MeV) (MeV) 



Single Breit-Wigner fit 



(1,1) 


766 ± 5 


258 ± 19 




1, 


.yo 


1 


U. 


.U / 


U. 


.4oU 


(1,2) 


769 ± 12 


263 ± 45 




o 

/, 


Of! 


1 


U. 


1 7 
.1 1 


U. 


/lOS 


(2,1) 


766 ± 10 


255 ± 37 




z. 


1 n 


1 

± 


U. 


1 K 


U. 


OAf\ 
.Z4U 


(2,2) 


768 ± 12 


303 ± 49 




2, 


.48 


± 


0. 


.16 


0. 


.816 


Teit-Wigner fit with 


constant background 
















(1,1) 


767±5 


210 ± 43 


0.19 ± 0.17 


1, 


.69 


± 


0. 


.23 


0. 


.365 


(1,2) 


768 ± 12 


209 ± 99 


0.20 ± 40 


1. 


.92 


± 


0. 


.61 


0. 


.470 


(2,1) 


766 ± 9 


208 ± 82 


0.17 ± 0.32 


1, 


.91 


± 


0. 


.47 


0. 


.218 


(2,2) 


765 ± 10 


188 ± 76 


0.41 ± 0.34 


1, 


.85 


± 


0. 


.42 


0. 


.700 



Table 7. Results of the fits to the four solutions of the S'-wave intensity measured 
in n~p 'K~'K^n at 17.2 GeV/c using Monte Carlo method for amplitude analysis. 
The fits are made with Breit-Wigner parametrization (6.3) and (6.4) with the Pisut- 
Roos shape factor. 
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Is B Ns xVdof 

Solution (MeV) (MeV) 



Single Breit-Wigner fit 

(1.1) 760 ± 8 269 ± 29 - 2.00 ± 0.13 0.414 

(2.2) 786 ± 21 408 ± 90 - 2.24 ± 0.16 1.140 
Breit-Wigner fit with constant background 

(1.1) 761±8 227 ± 68 0.12 ± 0.19 1.84 ± 0.28 0.394 

(2.2) 780 ± 13 187 ± 77 0.63 ± 0.31 1.51 ± 0.31 0.864 



Table 8. Results of the fits to two of four solutions of the 5'-wave intensity measured 
in n~p n~'K~^n at 17.2 GeV/c using minimization method for amplitude 
analysis. The fits are made with Breit-Wigner parametrization (6.3) and (6.4) with 
the Pisut-Roos shape factor. 
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Is B Ns xVdof 

Solution (MeV) (MeV) 



Single Breit-Wigner fit 



(1,1) 


723 ± 22 


282 ± 68 




U, 


.oo 


1 

± 


A 

u. 


1 n 
.lU 


U. 


QQQ 
.OOO 


(1,2) 


696 ± 36 


333 ± 128 




1, 




_i_ 
± 


n 
U. 




n 
U. 


.116 


(2,1) 


740 ± 32 


296 ± 116 




1, 


no 
.Uz 


1 

± 


U. 


on 

.zy 


U. 


.ZU4 


(2,2) 


714 ± 27 


362 ± 102 




1, 


.52 


± 


0. 


,30 


0. 


.194 


Teit-Wigner fit with 


constant background 
















(1,1) 


746±16 


145 ± 69 


0.18 ± 0.10 


0, 


.59 


± 


0. 


.18 


0. 


.712 


(1,2) 


706 ± 39 


262 ± 24 


0.13 ± 0.39 


1, 


.05 


± 


0. 


.44 


0. 


.114 


(2,1) 


745 ± 30 


165 ± 112 


0.23 ± 0.21 


0, 


.97 


± 


0. 


.42 


0. 


.094 


(2,2) 


724 ± 25 


211 ± 117 


0.25 ± 0.20 


1, 


.41 


± 


0. 


.42 


0. 


.124 



Table 9. Results of the fits to the four solutions of the S'-wave intensity measured 
in 7T~^n n~^'K~p at 5.98 GeV/c using Monte Carlo method for amplitude analysis. 
The fits are made with Breit-Wigner parametrization (6.3) and (6.4) with the Pisut- 
Roos shape factor. 
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Is B Ns xVdof 

Solution (MeV) (MeV) 



Single Breit-Wigner fit 



(1,1) 


778 ± 10 


158 ± 21 




1 1 o 

i.iy 


1 


U. 


1 o 


z. 


1 

.loo 


(1,2) 


749 ± 31 


353 ± 88 




1 /1 7 
1.4 ( 


it 


U. 


.AO 


U. 


A Qn 


(2,1) 


752 ± 20 


237 ± 52 




1 c;n 


1 

± 


U. 


.Z 1 


U. 


.844 


(2,2) 


749 ± 19 


309 ± 63 




1.92 


± 


0. 


,24 


0. 


.632 


Teit-Wigner fit with 


constant background 














(1,1) 


782±9 


117 ± 26 


0.08 ± 0.03 


1.13 


± 


0. 


.16 


1. 


.024 


(1,2) 


770 ± 24 


202 ± 74 


0.09 ± 0.04 


1.52 


± 


0. 


.32 


0. 


.080 


(2,1) 


763 ± 18 


153 ± 55 


0.12 ± 0.05 


1.46 


± 


0. 


.39 


0. 


.236 


(2,2) 


756 ± 15 


200 ± 59 


0.11 ± 0.05 


1.93 


± 


0. 


.33 


0. 


.212 



Table 10. Results of the fits to the four solutions of the S-wave intensity measured 
in 7T~^n n'^n~p at 11.85 GeV/c using Monte Carlo method for amplitude analysis. 
The fits are made with Breit-Wigner parametrization (6.3) and (6.4) with the Pisut- 
Roos shape factor. 
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Figure Captions. 



Fig. 1. Mass dependence of unnormalized amplitudes I^SpE and jS'pE measured 
in TT-pt TT-TT+n at 17.2 GeV/c at -t = 0.005 - 0.20 (GeV/c)^ using the Monte 
Carlo method for amplitude analysis (Ref. 23). Both solutions for the amplitude 
resonate at 750 MeV while the amplitude \S\'^'E is nonresonating in both 
solutions. 

Fig. 2. Mass dependence of unnormalized amplitudes \S\'^T, and measured 
in TT-p^ TT-TT+n at 17.2 GeV/c at -t = 0.005 - 0.20 (GeV/c)^ using the 
minimization method for amplitude analysis. Based on Fig. 10 of Ref. 13 and 
Fig. VI-21 of Ref. 12. Both solutions for the amplitude |5'pE resonate at 750 MeV 
while the amplitude |S'pE is nonresonating in both solutions. The analysis used 
the same data as in Fig. 1 (20 MeV mass bins). 

Fig. 3. The fits to amplitude jS'pS using the single Breit-Wigner parametrization 
(4.1) with Pisut-Roos shape factor (4.2a). The fitted parameters are given in Table 
1. 

Fig. 4. The fits to amplitude \S\'^T, using the single Breit-Wigner parametrization 
(4.1) with phenomenological shape factor F = 1. The fitted parameters are given 
in Table 1. 

Fig. 5. The fits to amplitude using the Breit-Wigner parametrization (4.6) 

with constant incoherent background and with phenomenological shape factors F = 
1. The fitted parameters are given in Table 2. 

Fig. 6. The fits to amplitude \S\^Ti using the Breit-Wigner parametrization (4.15) 
with constant coherent background and with phenomenological shape factor F = 1. 
The fitted parameters are given in Table 3. 

Fig. 7. The fits to amplitude jS'pS using the Breit-Wigner parametrization (4.19) 
with t-averaged constant coherent background and with phenomenological shape 
factor F — 1. The fitted parameters are given in Table 4. 
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Fig. 8. Mass dependence of unnormalized amplitudes |5'pS and measured in 

TT-p TT-TT+n at 17.2 GeV/c at -t = 0.005 - 0.20 (GeV/c)^ in 40 MeV mass bins 
from 600 to 1520 MeV. Based on Figs. 2 and 6 from Ref. 14. The amplitude 
resonates at 750 in Solution 1 and at 800 MeV in Solution 2 while the amplitude 
for |5'pS is nonresonating in both solutions in this mass range. 

Fig. 9. The fits to amplitude \S\^Ti using the Breit-Wigner parametrization (5.4) 
below 1120 MeV and a single Breit-Wigner formula with incoherent background 
above 1120 MeV. The phenomenological shape factor F = 1. The fitted parameters 
are given in Tables 5 and 6. 

Fig. 10. Four solutions for the S'-wave intensity Is measured in the reaction 
TT-p^ TT-TT+n at 17.2 GeV/c and -t = 0.005 - 0.2 (GeV/c)^ using Monte 
Carlo method for amplitude analysis (Ref. 23). The solid curves are fits to sin- 
gle Breit-Wigner parametrization (6.3). The dashed curves are fits to Breit-Wigner 
parametrization (6.4) with incoherent background. The fitted parameters are given 
in Table 7. 

Fig. 11. Two of the four solutions for the S'-wave intensity Is measured in the 
TT-p^ TT-TT+n at 17.2 GeV/c and -t = 0.005-0.2 (FeV/cf using minimization 
method for amplitude analysis. The data are based on Fig. 14a of Ref. 13 and Fig. 12 
of Ref. 11. The solid and dashed curves are Breit-Wigner fits as in Fig. 10. The 
fitted parameters are given in Table 8. 

Fig. 12. Four solutions for the S-wscve intensity Is measured in TT+n-^ — > TT+TT-p at 
5.98 GeV/c and —t = 0.2 — 0.4 (GeV/c)^ using Monte Carlo method for amplitude 
analysis (Ref. 23). The solid and dashed curves are Breit-Wigner fits as in Fig. 10. 
The fitted parameters are given in Table 9. 

Fig. 13. Four solutions for the S-wave intensity Is measured in TT+n^ — > TT+TT-p at 
11.85 GeV/c and —t = 0.2 — 0.4 (GeV/c)^ using Monte Carlo method for amplitude 
analysis (Ref. 23). The solid and dashed curves are Breit-Wigner fits as in Fig. 10. 
The fitted parameters are given in Table 10. 

Fig. 14. The t-evolution of mass dependence of moduli squared of t-channel nor- 
malized transversity amplitudes |Lp, |Lp, and |^7p in 7r+n^ — > TT+Tr~p at 
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5.98 GeV/c together with results for ir'p^ tt-tt+ti at 17.2 GeV/c and t = 0.068 
(GeV/c)2. 

Fig. 15. The ratio of amphtudes with recoil nucleon trans versity "down" and "up" 
with dimeson helicity A = 0. The deviation from unity shows the strength of Ai- 
exchange amplitudes. Based on Fig. 6 of Ref. 14. In our notation, gs = S, hs = S, 
gp = L, hp = L. 

Fig. 16. Test of predictions p^^ + pgo + ^Pn = —'^{Poo ~ Pii) = +2Pi-i due to 
vanishing ^i-exchange in 7r~p-|- — > tt^tt+tt- at 17.2 GeV/c and —t = 0.005 — 0.2 
(GeV/c)2. 

Fig. 17. Test of predictions Rep^Q = Repfg = Repg^ = due to vanishing Ai- 
exchange in ir-p-f tt-tt+ti at 17.2 GeV/c and -t = 0.005 - 0.2 (GeV/c)^. 

Fig. 18. S-wsive intensity normalized to 1 at maximum value. The data correspond 
to solutions Is{l, 1) and Is{2, 2) at 17.2 (GeV/c) from Ref. 23. The smooth curves 
are predictions of phase shift analysis for tt+tt~ — > tt+tt~ from Ref. 27. The dashed 
curve is the accepted solution Down, the dot-dashed curve is the rejected solution 
Up. 
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